DUDLEY  KNOX  LIBRARY 

NAVAt  POSTGRADUATE  SCHOOL  l 

MONTEREY  CA  93943-5101  us 


,  t. 


DEVELOPMENT  OF  A  COMPUTER  PROGRAM  FOR 
DETERMINATION  OF  SEARCH  AREAS  FOR  SEARCH  AND  RESCUE 
OPERATIONS  OF  THE  UNITED  STATES  COAST  GUARD 


by 

David  F.  Cunningham 

Lieutenant,  United  States  Coast  Guard 

and 

Donald  L.  Frantz 

Lieutenant  Commander,   United  States  Coast  Guard 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 
MANAGEMENT 

United  States  Naval  Postgraduate  School 
Monterey,  California 

19  6  5 


VPS  IWcW< 
WfcC 


CX>£ 


Libra 

U.  S.  Naval  Pi 
Mor' 


DEVELOPMENT  OF  A  COMPUTER  PROGRAM  FOR 
DETERMINATION  OF  SEARCH  AREAS  FOR  SEARCH  AND  RESCUE 
OPERATIONS  OF  THE  UNITED  STATES  COAST  GUARD 

by 
David  F.  Cunningham 
and 
Donald  L.  Frantz 
This  work  is  accepted  as  fulfilling 
the  research  paper  requirements  for  the  degree  of 
MASTER  OF  SCIENCE 
IN 
MANAGEMENT 
from  the 
United  States  Naval  Postgraduate  School 


ABSTRACT 

The  increase  in  trans  -oceanic  travel  and  maritime  operations  have 
created  a  need  for  new  and  improved  search  and  rescue  techniques  and 
facilities.    This  study  briefly  describes  the  various  phases  of  search 
and  rescue  operations  and  reviews  the  doctrine  currently  used  by  this 
Nation  as  set  forth  in  the  National  Search  and  Rescue  Manual. 

Attention  is  then  directed  to  developing  a  computer  program  which 
utilizes  current  doctrine  in  the  solution  of  one  portion  of  the  search  and 
rescue  problem.    The  program  determines  the  most  probable  location  of 
survivors  (datum)  and  the  associated  search  areas  at  any  specified 
time  for  an  off-shore  distress  case  with  an  estimated  reported  distress 
incident  position.    The  resulting  program/  written  in  Fortran  60,  is 
included  with  sample  test  data. 

Finally,  recent  developments  in  computer-generated  informations 
systems  which  may  offer  improvements  in  the  existing  doctrine  are 
discussed. 
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CHAPTER  I 
STATEMENT  OF  THE  PROBLEM 

Man's  continuing  struggle  against  the  elements  at  sea  is  noted 
throughout  history,     In  the  Book  of  Psalms  it  is  recorded: 

They  that  go  down  to  the  sea  in  ships,  that  do  business 
in  great  waters; 

These  see  the  works  of  the  LORD,   and  his  wonders  in 
the  deep, 

For  he  commandeth,   and  raiseth  the  storm  wind,  which 
lifteth  up  the  waves  thereof. 

They  mount  up  to  the  heaven,  they  go  down  again  to 
the  depths:  their  soul  is  melted  because  of  trouble. 

They  reel  to  and  fro,  and  stagger  like  a  drunken  man, 
and  are  at  their  wit's  end,  1 

Each  jiear  literally  thousands  of  people  find  themselves  in  distress 

while  located  upon  the  high  seas  or  upon  waters  over  which  the  United 

States  has  jurisdiction.    Fortunately  for  most  of  them,  their  rescue  is 

only  a  matter  of  waiting  a  relatively  short  time.     This  life-saving 

assistance  usually  appears  in  the  form  of  one  or  more  units  operated  by 

the  United  States  Coast  Guard,    Specifically  charged  by  Title  14, 

United  States  Code.  *  the  United  States  Coast  Guard  carries  the 

responsibility  of  establishing^  operating  and  maintaining  these  facilities 

^•The  Holy  Bible  (Kino  James  Version;  Cleveland:    The  World  Publishing 
Company,  n.d.)#   Psalms  107:23-27, 

United  States  Code,   Title  14  (Coast  Guard  Chapter  1,  Section  2, 
Vol,  III,    (1958  edition);  Washington:  Government  Printing  Office,    1959), 
p,   2286. 


necessary  to  provide  such  assistance  as  may  be  required  by  any  person 
or  persons,   vessels  or  aircraft  in  distress  upon  these  waters. 

I.      THEORETICAL  BACKGROUND 

The  task  of  successfully  conducting  search  and  rescue  operations 
at  sea  has  been  a  challenge  ever  since  man  first  ventured  out  upon  the 
hydrosphere.     Modern  technology  has  contributed  greatly  to  providing 
the  answer  to  this  challenge.     Nevertheless,   each  individual  case 
presents  the  would-be-rescuers  with  many  complex  problems  faught  with 
uncertainties  and  the  capaciousness  of  Nature. 

To  increase  man's  effectiveness  in  rendering  assistance  to  his 
fellow  man  on  the  high  seas,   the  martime  nations  of  the  world  have 
combined  their  efforts  in  developing  and  improving  search  and  rescue 
techniques,   establishing  necessary  communication  systems,  and 
providing  for  mutual  assistance  and  cooperation  in  large  scale  operations, 
The  National  Search  and  Rescue  Manual,  CG-308,   developed  by  the 
United  States  Coast  Guard,   serves  as  this  Nation's  leading  guide  in 
this  area. 3 

° United  States  Coast  Guard,  National  Search  and  Rescue  Manual, 
CG-308  (Washington:  Government  Printing  Office,    1959,  Amended  1964). 


Generally,   search  and  rescue  operations  can  be  divided  into  four  main 
phases.     The  first  of  these  is  the  "Alert"  when  a  distress  or  potential 
distress  situation  is  suspected.     This  would  occur,   for  example,  when 
an  aircraft  or  boat  is  overdue  or  has  failed  to  meet  routine  scheduled 
communications  checks.     During  this  phase  efforts  are  directed  to 
determine  if,   in  fact,  there  is  an  overdue  craft,  and  whether  or  not 
communications  can  be  established  with  the  craft  in  question. 

The  second,  the  "Distress"  phase,  follows  when  the  "Alert"  fails 
to  establish  the  fact  that  a  distress  does  not  exist.     Here  one  finds 
efforts  being  directed  to  determine  the  most  likely  location  of  the 
distressed  unit  with  an  associated  area  for  searching  and  the  assignment 
and  dispatching  of  search  units  to  the  scene  of  the  search. 

The  third  is  the  "Search"  phase  wherein  those  search  units  previously 
designated  conduct  searches  of  the  area  determined  in  the  second  phase. 

The  fourth  and  last  phase  is  the  "Rescue"  which  hopefully  follows 
a  successful  "Search.  " 

It  is  evident  when  a  distress  incident  does  exist  that  no  amount  of 
searching  will  lead  to  rescue  if  this  searching  has  been  conducted  in  an 
area  which  does  not  contain  the  unit  in  distress.    Therefore,  it  is 
essential  that  the  second  phase  produce  sound  information  and  directions 
for  the  search  units.     This  requires  the  careful  analysis  and  evaluation  of 


many  complex  and  interrelated  factors.     Some  of  these  factors  are 

nature  of  the  distress,   reliability  of  navigational  information  available 

on. the  probable  position  of  the  distressed  unit,  the  effects  of  ocean  currents, 

local  wind-driven  currents,   sea  and  swell  conditions,   life  raft  drift  or 

leeway,   possible  effect  of  parachute  drift  if  such  is  present  in  the 

distress  incident  under  consideration,   estimated  time  of  survival  under 

conditions  existing  at  the  scene,  time  involved  for  the  search  units  to 

reach  the  scene  with  the  associated  change  in  site  location  caused  by 

the  factors  stated  earlier,  the  navigational  errors  of  the  search  units, 

and  the  various  probabilities  associated  with  detection  in  different  size 

areas. 

At  the  present  time  the  evaluation  and  analysis  of  these  factors  are 
carried  out  by  human  decision-makers  with  the  ai$  of  excellent  guide- 
lines,  charts,   graphs  and  tables.    The  final  decision  is,   therefore,  a 
human  one  which  requires  time  to  develop,   possesses  inherent  chances 
for  error  in  the  calculations,  and,   by  its  very  nature,   must  be  subject 
to  the  human  bias  and  judgement  of  the  decision-maker  which  may  or 
may  not  adversely  effect  the  decision. 

II.      STATEMENT  OF  SPECIAL  PURPOSE 

The  advent  of  the  computer  has  greatly  facilitated  scientific  analysis 
and  routine  decision-making.     To  date  the  computer  has  not  been  applied 


to  the  search  and  rescue  problem  although  the  idea  has  been  expressed 
and  has  been  the  subject  of  some  feasibility  studies.     The  application 
of  computers  in  this  field  should  greatly  facilitate  the  analysis  of  the 
many  factors  involved  and  also  reduce  the  time  required  to  reach  a 
decision.     In  view  of  the  fact  that  human  lives  are  at  stake  in  these 
operations,   it  is  mandatory  that  the  best  possible  decisions  be  made 
with  the  information  and  alternatives  available,  and  these  decisions 
must  be  made  in  the  shortest  possible  time.     The  following  problem  has 
been  selected  with  these  thoughts  in  mind. 

III.       THE  PROBLEM 

It  is  the  purpose  of  this  study  (1)  to  completely  develop  a  computer 
flow  diagram  and  associated  computer  program  using  the  Control  Data 
Corporation  compiler  Fortran  60  language  (equivalent  to  Fortran  II)  for 
that  portion  of  the  "Distress"  phase  of  search  and  rescue  (SAR)  operations 
relating  to  the  determination  of  a  datum  point  and  associated  search  area 
based  on  the  doctrine  outlined  in  the  National  Search  and  Rescue  Manual, 
Chapter  Six;  and  (2)  to  investigate  the  possible  development  of  those 
programs  encountered  during  the  research  process  which  appear  to  offer 
improved  techniques  for  analyzing  and  evaluating  various  portions  of  the 
problems  encountered  in  the  "Distress"  phase.     These  potentially  improved 
techniques  will  be  submitted  to  the  United  States  Coast  Guard  for 
evaluation. 


IV.       DEFINITIONS  OF  KEY  TERMS 

The  major  key  term  under  consideration  is  the  "Distress"  phase.    As 
used  in  Chapter  Six  of  the  National  Searph  and  Rescue  Manual  this  refers 
to  that  phase  of  search  and  rescue  operations  wherein  the  estimated 
location  of  the  distress  incident  and  an  associated  surrounding  area 
called  the  "search  area"  are  determined  and  designated  for  any  specified 
time.     It  also  includes  the  assignment  and  dispatching  of  search  units  to 
the  scene  of  the  distress  incident.     However,   for  the  purpose  of  this 
study  no  consideration  will  be  given  to  the  assignment  and  dispatching 
of  search  units. 

The  following  terms  directly  associated  with  the  "Distress"  phase 
will  be  used: 

Distressed  unit.    The  distressed  unit  is  that  craft  whose  survivors 
require  rescue  or  assistance. 

The  object  of  the  search  may  consist  of  a  disabled  vessel  of  any 
type,   a  downed  aircraft,  a  life  raft  or  life  boat,   or  any  form  of  survival 
equipment  which  will  support  a  survivor  or  survivors  on  the  water's 
surface,  or,   in  some  cases,   the  bodies  of  those  who  did  not  survive  but 
which  may  be  floating. 

Rescue  Coordination  Center.    The  Rescue  Coordination  Center  (RCC) 
is  that  operational  office  responsible  for  the  initiation,   coordination  and 


overall  control  of  search  and  rescue  operations  for  the  specific  areas 
under  its  jurisdiction.     For  Pacific  operations  this  is  controlled  by 
Commander,  Western  Area,  and  for  Atlantic  operations  it  is  controlled  by 
Commander,   Eastern  Area.     The  specific  individual  responsible  for  the 

various  decisions  required  on  any  given  distress  case  is  referred  to  as 
the  RCC  Controller. 

Qn-Scene  Commander.    This  is  the  commander  of  a  search  unit  at 
or  near  the  scene  of  the  search  who  has  been  designated  by  the  RCC 
Controller  as  the  commander  responsible  for  controlling  and  coordinating 
the  efforts  of  the  various  units  participating  in  the  search  while  they  are 
operating  at  the  scene. 

Search  unit.  Any  unit  which  is  actively  participating  in  the  search 
while  under  some  form  of  control  of  the  Rescue  Coordination  Center  (RCC) 
or  On-Scene  Commander  (OSC)  is  designated  as  a  search  unit.     This  may 
be  an  aircraft  or  surface  vessel,   commerical  or  military. 

Distress  incident  position.    This  is  the  geographical  position  of  the 
distressed  unit  at  the  time  that  the  distress  ocurred.     The  accuracy  with 

which  this  position  is  determined  will  depend  upon  the  specific  nature  of 
the  distress  and  the  time  and  facilities  available  for  determining  the 
position. 
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Datum.    The  datum  point  is  that  geographical  position  which  represents 
the  most  likely  location  of  the  survivors  at  any  specified  time.     Generally, 
the  initial  datum  will  correspond  with  the  "distress  incident  position.  " 
At  any  time  after  the  distress  has  occurred,   the  datum  point  will  have 
moved  in  response  to  the  forces  of  Nature  and  the  efforts  of  the  survivors. 
The  datum  point  is  considered  to  be  in  the  center  of  any  designated  search 
area0     It  serves  as  a  reference  point  for  determining  the  limits  of  a 
designated  search  area. 

Search  area.    That  area  in  which  the  survivors  are  believed  to  be 
located  and  which  can  be  covered  with  a  reasonable  degree  of  thoroughness 
by  the  search  units  is  defined  as  the  search  area.     It  is  bounded  by  limits 
established  by  the  search  radius  which  extends  outward  from  the  datum 
point. 

Search  radius.    The  search  radius,    (R),   is  the  distance  (in  nautical 
miles)  from  the  datum  point  to  the  edge  of  the  search  area.     It  is  a 
function  of  the  number  of  the  search,   for  example--the  second  search, 
and  the  total  probable  error  of  position.     The  number  of  the  search  determines 
the  safety  factor  to  be  applied.     The  total  probable  error  of  position  is  a 
function  of  several  other  factors  which  are  defined  below.     In  general,   if 
the  search  radius  is  R,   the  safety  factor  is  SF,  and  total  probable  error 
of  position  is  c,   then 

R  =  SF  x  c 
8 


Total  probable  error  of  position.    The  total  probable  error  of  position, 
(c),   is  that  distance  (in  nautical  miles)  which  represents  the  best 
estimate  of  position  error  when  the  following  factors  are  taken  into 
consideration: 

a.  Initial  position  error  of  the  distressed  unit,   (X). 

b.  Navigational  error  of  the  search  unit(s),    (Y). 

c.  Life  raft  drift  error,    (de). 

Total  probable  error  of  position  is  a  scalar  quantity  which  is  applied  to 
vectors  extending  in  all  directions  from  the  datum  point.     It  is  represented 
by  the  equation: 

c  =  Vde2  +  X2  +  Y2 
A  review  of  the  development  of  this  equation  is  contained  in  Chapter  II. 

Initial  position  error,   (X).    This  is  an  estimate  of  the  distance  (in 
nautical  miles)  by  which  the  distress  incident  position  may  be  in  error. 
It  is  a  function  of  the  method  of  navigation  used  to  fix  the  position  of  the 
distressed  craft. 

Navigational  error  of  the  search  unit,   (Y).    The  navigation  used  in 
guiding  the  search  units  to  the  distress  incident  position  or  datum  point 
and  throughout  the  conduction  of  the  search  will  be  subject  to  varying 
degrees  of  error  depending  upon  the  method  of  navigation  used.    This 
error  (in  nautical  miles)  must  be  considered  when  determining  the  total 
probable  error  of  position,   c,  and  the  resulting  search  radius. 


Life  raft  drift  error,  fcU).  Once  survivors  are  free  of  the  distressed  unit 
and  floating  upon  the  surface,   they  will  be  subjected  to  both  surface  currents 
and  surface  winds.     The  interaction  of  these  forces,   the  nature  of  the 
survival  craft--that  is,  whether  it  is  a  life  raft,   life  boat  or  life  jacket — 
and  the  time  span  over  which  they  react,  all  contribute  to  a  resultant 
drift  vector,    D.     This  drift  vector  may  have  maximum  and  minimum  values 
depending  upon  the  variations  in  the  forces  of  Nature  during  the  period  of 
interaction  and  the  availability  of  accurate  information  about  the  weather 
conditions  at  the  scene  of  the  distress  incident.     This  variation  introduces 
the  concept  of  life  raft  drift  error,    (de),  which  will  contribute  to  the  total 
probable  error  of  position. 

Drift.    The  movement  of  the  survivors  (datum)  upon  the  surface  is 
represented  by  the  vector,    D.     This  is  the  vector  which  is  the  sum  of  the 
average  sea  current  vector,   the  local  wind  current  vector,   and  the  leeway 
vector.    Although  the  average  sea  current  is  assumed  to  have  a  constant 
value,    the  latter  two  vectors  may  vary  thus  contributing  to  maximum  and 
minimum  values  for  drift,    D. 

Drift  (max).   Drift  (max),    Dmax,   is  that  vector  which  represents  the 
greatest  displacement  possible  for  a  given  datum  during  a  specified 
period  of  time,  At. 

Drift  tmin).    Drift  (min),   Dm^n,   is  that  vector  which  represents  the 

minimum  displacement  possible  for  a  given  datum  during  a  specified 

period  of  time,  At. 
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FIGURE    1 

COMPONENTS  OF  DRIFT,    D,   WHEN  FLUCTUATIONS  ARE 
SMALL  OR  EXISTING  CONDITIONS  ARE  KNOWN 
WITH  A  HIGH  DEGREE  OF  ACCURACY 
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FIGURE  2 

A  MINIMAX  PLOT  OF  DRIFT  COMPONENTS  WHEN  FLUCTUATIONS 
ARE  L^RGE  OR  UNCERTAINTY  EXISTS  AS  TO  THE 
EXACT  CONDITIONS  AT  THE  SCENE 
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Average  sea  current.    The  average  sea  current  is  that  vector  described 
by  the  set  and  drift  of  the  prevailing  ocean  currents  for  the  area  of  the 
distress  at  the  time  of  the  incident.     It  is  determined  by  extracting  the 
applicable  values  from  the  appropriate  charts  in  the  Atlas  of  Surface 
Currents 4  for  the  area  to  be  searched.     The  average  sea  current  vector  is 
only  considered  if  the  Atlas  of  Surface  Currents  indicates  that  the  reported 
current  has  a  steadiness  equal  to  or  greater  than  35%  of  the  time. 

Local  wind  current.     The  direction  and  magnitude  of  the  average  sea 
current  may  be  altered  by  the  presence  of  local  surface  winds,     The 
magnitude  of  the  current  generated  by  the  local  surface  winds  is  a 
function  of  the  force  of  the  winds,  their  steadiness  both  in  magnitude  and 
direction,  the  size  of  the  area  over  which  they  have  been  blowing  (fetch), 
the  length  of  time  that  they  have  been  blowing,  and  the  influence  of  the 
Coriolis  force. 


Atlas  of  Surface  Currents  (Washington:  U.   So   Navy  Hydrographic 
Office,    1953).     This  Atlas  consists  of  a  series  of  different  atlases  for 
various  ocean  areas  of  the  world.     For  any  given  area  the  charts  are 
furnished  for  each  month  of  the  year.     These  charts  are  divided  into  grids 
for  each  degree  of  latitude  and  longitude.    An  arrow  and  numbers  placed 
in  each  one-degree  quadrangle  shows  the  mean  direction  of  force  of  the 
surface  current  in  that  quadrangle  for  that  particular  month  under  average 
normal  conditions  and  the  number  of  observations  reported  for  that 
quadrangle.    A  current  rose  shows  the  frequency  distribution  (steadiness) 
and  average  drift  rate  for  the  eight  cardinal  and  intercardinal  points  for 
larger  quadrangles. 
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Leeway.   In  addition  to  generating  local  wind  currents,  the  local 
surface  winds  also  act  directly  upon  that  portion  of  the  survivors  or  their 
survival  craft  which  is  above  the  surface.     The  resulting  displacement  is 
called  "Leeway.  " 

Parachute  drift.     If  the  distress  incident  involves  an  aircraft  distress 
which  has  required  that  the  survivors  bail  out,   consideration  must  be  given 
to  their  free  fall  and  parachute  drift.     The  distress  incident  position  will 
be  the  position  of  the  aircraft  at  the  time  of  bail-out,  and  the  initial 
datum  will  be  that  position  where  the  parachute  would  have  landed  after 
considering  the  altitude  at  which  the  aircraft  was  flying  and  the  prevailing 
winds  at  different  altitudes  from  the  point  of  bail -out  to  the  surface. 

IV.      ASSUMPTIONS 

In  approaching  the  first  part  of  the  problem  it  has  been  assumed  that 
the  procedures,   graphs,   tables  and  guidelines  outlined  in  the  National 
Search  and  Rescue  Manual  are  based  on  the  best  information  and  theories 

available  for  the  Practical  methods  presently  used  by  the  United  States 
Coast  Guard  in  search  and  rescue  operations.     No  further  efforts  have 
been  directed  toward  verifying  or  disproving  the  existing  doctrine* 
Instead,   efforts  have  been  concentrated  on  the  development  of  an 
operational  computer  program  which  will  facilitate  the  application  of  the 
existing  doctrine. 
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In  developing  the  computer  program  it  was  assumed  that  the  various 
position  and  navigational  errors  could  be  approximated  by  normal 
distribution  curves.     Furthermore,  where  a  vector  produced  two  possible 
positions  (see  Figure  2,   page  11),   that  is  a  maximum  and  a  minimum 
position,   it  was  assumed  that  the  most  likely  position  would  be  at  a  point 
midway  between  the  two  positions  lying  on  a  line  connecting  the  two  points. 

V.      LIMITATIONS 

Certain  self-imposed  limitations  were  made  in  considering  the 
problem,     First,  the  problem  was  limited  to  the  type  where  the  time  that 
the  distress  occurred  was  known,   the  nature  of  the  distress  and  the  type 
of  units  involved  were  known,  and  an  estimated  distress  incident  position 
was  known.     Furthermore,   the  estimated  distress  incident  position  was 
placed  sufficiently  far  offshore  to  eliminate  the  effect  of  tidal  currents, 
rotary  currents  resulting  from  coastal  configuration,  and  currents  caused 
by  rivers  flowing  into  the  sea<, 

Additional  limitations  precluded  the  development  of  the  parachute 
drift  program.     Recent  research  into  this  area  has  resulted  in  the  development 
of  a  computer  program  to  handle  this  problem.    Where  an  aircraft  is 
involved  as  the  distressed  unit,   the  program  is  limited  to  the  type  of 
distress  which  involves  ditching  the  aircraft  into  the  ocean, 

No  attempt  was  made  to  develop  a  program  involving  a  small  craft 
as  the  distressed  unit  since  such  a  craft  is  unlikely  to  venture  far  enough 
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away  from  the  coast  to  fall  within  the  parameters  adopted  in  limiting  the 
distress  incident  position. 

Testing  of  the  model  was  done  by  generation  of  hypothetical  cases, 
manually  computing  a  solution  using  the  techniques  and  procedures 
outlined  in  the  National  Search  and  Rescue  Manual,  and  then  comparing 
these  results  with  those  produced  on  the  computer.     Historical  data  of 
actual  distress  cases  were  not  avajyLable  at  the  time  that  this  project 
was  completed. 

The  program  was  developed  for  operations  in  the  Atlantic  Ocean, 
north  of  the  equator  only.    With  minor  modifications,   the  program  can  be 
adapted  for  operations  south  of  the  equator  and  in  the  Pacific  Ocean,   both 
east  and  west  of  the  International  Date  Line. 

One  final  limitation  on  the  study  relates  to  the  method  used  in 
measuring  distances.    All  distances  are  treated  as  rhumb  line  distances 
based  on  Mercator  Sailings  rather  than  Great  Circle  distances.     Since 
the  magnitude  of  the  distances  involved  is  small,   the  error  resulting  in 
using  Mercator  distances  will  be  neglected. 

VI.      RESEARCH  SIGNIFICANCE 

The  development  of  an  operational  computer  program  to  assist  the 
decision-makers  responsible  for  coordinating  and  conducting  search  and 
rescue  operations  will  constitute  a  major  step  forward  in  improving  Coast 
Guard  operations.     By  eliminating  the  complex  and  tedious  manual 
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calculations  and  plotting  of  information,  it  will  free  the  Rescue 
Coordination  Center  Controller  and  his  team  to  concentrate  more  fully  on 
those  phases  of  the  operation  which  require  the  all  important  .tasks  of 
human  decision-making=-the  final  value  judgments. 

The  development  of  this  program  opens  new  horizons  for  the  use  of 
computers  in  expanding  search  and  rescue  operations  by  presenting  a 
workable  program  which  may  be  adapted  to  the  existing  Atlantic  Automated 
Merchant  Vessel  Reporting  System  (AMVER)  and  the  newly  developing 
Pacific  AMVER  system. 

The  most  important  significance  of  such  a  program  is  that  it  should 
lead  to  improved  effectiveness  in  that  mission  of  the  United  States  Coast 
Guard  which  is  directed  to  the  preservation  and  safety  of  the  lives  and 
property  of  all  those  who  travel  over  and  upon  the  seas. 
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CHAPTER  II 
REVIEW  OF  THE  LITERATURE 

Studies  directed  toward  the  problem  of  search  and  rescue  have  been 
limited  primarily  to  those  conducted  by  the  various  Armed  Forces  and 
those  specific  agencies  of  other  countries  that  have  been  responsible 
for  the  execution  of  this  type  of  work.     The  United  States  Coast  Guard 
has  been  the  leader  in  this  field  both  in  the  United  States  and  in 
international  efforts  to  develop  the  best  techniques  possible.    Accordingly, 
the  literature  dealing  with  this  problem  has  been  produced  by  the  combined 
efforts  of  the  various  services  and  agencies  and  certain  individuals 
directly  connected  with  these  operations. 

Since  the  problem  of  search  and  rescue  is  amenable  to  study  by  the 
various  phases  involved,   this  review  of  the  literature  will  be  limited  to 
that  phase  dealing  with  the  determination  of  the  best  estimated  position 
of  the  survivors  at  any  given  time  having  started  with  an  initial  estimated 
position,   and  then  the  determination  of  a  reasonable  area  for  searching 
around  this  estimated  position,     Furthermore,   this  review  will  be  focused 
on  two  specific  areas.     The  first  will  be  devoted  to  the  doctrine  prescribed 
in  the  National  Search  and  Rescue  Manual,   Chapter  Six,    "Determination 
of  Search  Areas.  "    It  is  for  a  portion  of  this  Chapter  that  the  computer 
program  has  been  developed. 
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The  second  area  deals  with  the  literature  which  sheds  new  light  on 
certain  aspects  of  the  problem  which  may  tend  to  improve  the  existing 
doctrine  as  currently  outlined  in  the  National  Search  and  Rescue  Manual. 
Therefore,   all  literature  reviewed  will  be  that  which  has  been  developed 
reflecting  the  latest  theories  and  ideas  since  the  development  of  the 
National  Search  and  Rescue  Manual  in  1961. 

I.       DOCTRINE  ON  THE  DETERMINATION  OF  THE  SEARCH  AREA 

ACCORDING  TO  THE  NATIONAL  SEARCH  AND  RESCUE  MANUAL 

Before  search  and  rescue  units  are  dispatched  to  the  scene  of  a 
distress  and  before  the  type  of  pattern  of  search  to  be  used  is  determined, 
an  estimated  position  of  the  survivors  must  be  established,   and  an  area 
sufficiently  large  to  insure  inclusion  of  the  survivors  must  be  selected. 
This  search  area  is  a  function  of  the  various  parameters  of  the  incident 
and  the  associated  position  errors.     There  are  two  possible  situations 
confronting  the  decision-maker.     First,   the  initial  distress   position  is 
known  with  varying  degrees  of  accuracy;  and  the  second,   the  initial 
position  of  distress  is  completely  unknown.     Only  the  first  case  will  be 
covered  in  this  paper. 

There  are  many  factors  associated  with  determining  the  initial  distress 
position  and  its  related  position  errors.     The  following  are  some  of  the 
questions  which  must  be  answered  by  the  search  and  rescue  (SAR) 
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coordinator  (RCC  Controller).    What  type  of  craft  is  in  distress?    What 
is  the  nature  of  the  distress?    What  type  of  navigational  information  was 
used  to  determine  the  initial  distress  position?    What  type  of  survival 
craft  is  likely  to  be  employed  by  the  survivors?    What  are  the  existing 
weather  conditions  at  the  scene  of  distress?    What  has  the  weather  been 
like  at  the  scene  for  the  preceding  twenty-four  hours?    How  long  will  it 
be  before  a  SAR  unit  can  reach  the  scene?    Where  is  the  most  likely 
position  (datum)  of  th a  survivors  when  the  SAR  units  reach  the  scene? 
What  are  the  various  drift  factors  which  must  be  considered  in  this 
particular  case?    What  is  the  expected  time  of  survival  for  the  survivors 
considering  sea  water  temperature  and  surface  weather  conditions  at  the 
scene? 

In  approaching  this  problem,   the  National  Search  and  Rescue  Manual 
first  divides  the  factors  affecting  the  position  into  three  broad  classes. 
The  first  of  these  is  called  the  "Initial  Error  in  Position" ^  and  is  a  function 
of  both  the  type  of  craft  involved  and  the  method  of  navigation  used.     If 
the  reported  initial  position  of  the  distress  incident  is  based  upon  an 
accurate  navigational  fix,   the  following  position  errors  are  assumed:" 


National  Search  and  Rescue  Manual,   op.  cit. ,   p,   6-2, 


6Ibid. ,   p.   6-7 
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Radius  of 
Type  of  distress  craft  Assumed  Position  Error 

Ship  Fix  Five  nautical  miles 

Aircraft  Fix  Ten  nautical  miles 

Small  Craft  Fix  Fifteen  nautical  miles 

If  the  initial  position  of  the  distress  incident  was  determined  by  monitor 

stations  on  radar  nets  or  on  military  or  commerical  radio  direction-finding 

nets,   these  various  nets  will  classify  their  fix.     These  classifications 

have  standard  specified  position  errors. 

If  a  fix  could  not  be  obtained  at  the  time  of  the  incident,   or  if  the 

reported  position  was  not  specifically  reported  as  a  fix,   it  is  assumed  that 

the  position  is  based  on  dead  reckoning  (DR).     In  this  case  the  following 

position  errors  are  assumed. 

Type  of  distress  craft  Radius  of  Assumed  Position  Error 

Five  percent  of  the  distance 
Ship  DR  traveled  since  the  last  fix, 

plus  the  error  of  that  fix. 

Ten  percent  of  the  distance 
Aircraft  DR  traveled  since  the  last  fix, 

plus  the  error  of  that  fix. 

Fifteen  percent  of  the  distance 
Small  Craft  DR  traveled  since  the  last  fix, 

plus  the  error  of  that  fix. 

The  second  factor  affecting  the  distress  position  is  that  of  parachute 

drift  where  the  particular  distress  incident  involves  an  aircraft  foom  which 

the  survivors  have  bailed  out.     The  position  reported  as  the  initial 


7  Ibid. 
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position  of  the  distress  incident  will  not  be  the  same  position  as  that  of 
the  survivors  when  they  land  in  the  water,,    A  table  is  provided  to  determine 
how  far  downwind  the  landing  will  be  from  the  initial  reported  position  of 
distress  or  bail-out  point.  The  entering  arguments  for  this  table  are 
altitude  of  parachute  opening  and  the  wind  velocity  in  knots  which 
represents  the  average  wind  between  the  winds  at  the  altitude  at  which  the 
parachute  opens  and  the  surface,  ° 

The  third  factor,   and  the  one  most  subject  to  question,   is  that  of 
"Survival  Craft  Drift.  "9    Once  the  survivors  are  adrift,   they  are  subjected 
to  the  various  elements  of  Nature.     In  the  open  sea  the  simplifying 
assumption  has  been  made  that  the  resultant  movement  of  the  survivors 
is  the  consequence  of  three  forces  acting  upon  them.     These  three  forces 
influence  the  survivors  motion  or  travel  independently.     They  are  (1)  the 
force  due  to  the  average  sea  current  which  is  a  function  of  the  major 
ocean  current  circulation  pattern;  (2)  the  force  due  to  local  wind  currents 
which  have  been  generated  by  the  local  winds  in  the  area  for  the  twenty- 
four  hours  preceding  the  distress  incident;  and  (3)  the  leeway  of  the 
survival  craft  as  it  responds  to  the  local  winds  at  and  following  the  time 
of  the  incident. 


8Ibid. ,    p.   6-3, 


Ibid.      />, 
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The  average  sea  current  is  determined  by  consulting  the  appropriate 
Atlas  of  Surface  Currents  for  the  area  and  month  involved.     If  the  average 
sea  current  is  indicated  as  having  a  steadiness  (frequency  of  reports)  of 
35%  or  greater,   it  is  considered  as  a  relevant  force.     Otherwise,  it  is 
disregarded. 

The  local  wind  current  is  treated  as  follows: 

a.  Obtain  an  estimate  of  the  local  wind  vector  for  the  preceding 
twenty-four  hour  period.     If  the  vector  has  not  varied,  assume  that  the 
local  wind  current  has  diverged  from  the  wind  vector  by  thirty  degrees  to 
the  right  for  latitudes  north  of  10°N,  thirty  degrees  to  the  left  for  latitudes 
south  of  10°S,   or  is  directly  downwind  for  latitudes  10°N  to  10°S.     This 
divergence  reflects  the  effect  of  the  Gbriolis  force. 

b.  If  the  wind  vector  has  varied  during  the  twenty-foliihour 
period,   determine  the  mean  wind  for  four  hour  periods  during  the  twenty- 
four  hours.     Then  use  the  mean  wind  to  determine  the  direction, of  the 
local  wind  current.     If  variations  have  been  more  than  45°  or  ten  knots, 
intervals  smaller  than  four  hours  should  be  used,  with  the  final  local  wind 
current  direction  becoming  the  sum  of  the  vectorsi forcall the  time  intervals. 

c.  The  magnitude  of  the  local  wind  current  is  taken  from  a  table 
with  an  entering  argument  of  the  local  wind  speed  in  knots. 

d.  The  assumption  is  made  that  for  the  first  six  to  twelve  hours  that 
a  wind  is  blowing,   the  local  wind  current  generated  will  be  downwind.     The 
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divergence  indicated  in  a.  above  applies  only  after  the  wind  has  been 
blowing  for  more  than  twelve  hours. 

The  component  of  leeway  is  equally  complex.     For  survivors  supported 
by  life  jackets  or  floating  free  on  the  surface,   the  sail  area  exposed  to  the 
winds  is  considered  to  be  negligible  as  compared  to  the  submerged  surface 
area  which  will  be  subjected  to  the  surface  currents.     For  those  survivors 
supported  by  life  rafts,   the  leeway  will  be  directly  downwind.     The 
magnitude  of  this  leeway  will  be  dependent  upon  whether  or  not  the 
survivors  have  streamed  a  drogue.     If  the  survivors  are  in  a  life  boat,  the 
magnitude  will  also  be  dependent  upon  the  absence  or  presence  of  a  drogue. 
In  addition,   the  drift  of  a  life  boat  may  vary  as  much  as  forty  degrees  on 
either  side  of  a  downwind  direction  depending  upon  the  configuration  and 
loading  of  the  life  boat. 

If  the  winds  have  been  fluctuating,   the  leeway  is  determined  for  four- 
hour  periods  and  then  added  vectorially  to  determine  the  final  component. 
The  National  Search  and  Rescue  Manual  provides  a  graph  for  determining 
the  magnitude  of  a  typical  life  raft  leeway  (the  type  is  not  indicated)  for 
conditions  with  and  without  a  drogue.  ^ 

Whenever  uncertainty  exists  as  to  the  length  of  time  adrift,  the 
magnitude  of  the  leeway,   or  the  accuracy  of  the  reported  conditions  of 
wind  and  current,   leeway  is  approached  using  minimax  plotting  techniques.  ■*■* 


10Ibid.  .    p,   6-4, 
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Ibid. ,    p„   6-5. 
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Under  these  conditions  the  RCC  Controller  must  determine  the  minimum 

and  the  maximum  probable  values  of  wind  and  current.     From  this  information 

a  plot  is  constructed  to  reflect  the  total  minimum  possible  drift  and  the 

total  maximum  possible  drift.     The  "most  probable  location  of  the  target 

is  the  position  midway  between  the  minimum  and  maximum  drift  positions.  "^ 

If,   however,   the  wind  conditions  have  been  fairly  steady,  an  average 
value  may  be  used.     This  will  lead  to  a  final  value  of  the  "drift  vector  which 
is  then  considered  to  be  an  average  drift  vector.     Both  of  these  concepts 
were  illustrated  in  Figures  1  and  2,Chapter  I,   page  11. 

The  status  of  the  problem  at  this  point  finds  the  RCC  Controller  with 
an  estimated  position  of  the  incident  based  on  information  received  in  the 
"Alert"  phase.     He  also  can  determine  the  most  probable  location  of  the 
survivors  after  applying  the  drift  factors.     The  next  step  is  to  evaluate 
the  total  probable  error  associated  with  the  latter  position. 

The  total  probable  error  of  position  is  obtained  by  considering  the 
initial  position  error,   the  navigational  error  of  the  search  unit(s),   and  the 
drift  error  associated  with  the  determination  of  the  drift  of  the  survivors. 
The  first  of  these,   initial  position  error,   is  determined  in  the  manner 
described  earlier  in  this  Chapter.     The  navigational  error  of  the  search 
unit(s)  is  handled  in  a  similar  manner  using  an  identical  table  for 
determining  the  value  of  the  error.     The  values  for  these  two  errors  are 


12  Ibid, 
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represented  by  "X"  and  "Y"  respectively. 

The  drift  error,   de,   can  be  determined  one  of  three  ways.     If  the 
total  drift  was  determined  as  an  average  drift,  the  drift  error  is  estimated 
to  be  equal  to  the  total  drift  divided  by  eight.  "    if  minimum  and  maximum 
drifts  have  been  computed,  a  table  is  provided  for  determining  the  drift 
error  using  the  minimum  and  maximum  drift  values  as  entering  arguments. 

The  third  method  of  determining  drift  error  is  used  when  the  maximum 
and  minimum  drifts  have  been  plotted.     Using  the  values  for  the  maximum 
and  minimum  drift,   the  respective  drift  errors  are  determined  by  dividing 
each  value  by  eight.     These  values  for  de  are  then  plotted  as  the  radii  of 
circles  centered  on  their  respective  drift  positions.     The  final  drift  error, 
referred  to  as  the  minimax  drift  error,   is  found  graphically  to  be  the 
radius  of  a  circle  centered  on  the  line  connecting  the  maximum  and 
minimum  drift  positions  and  which  is  tangent  to  the  outer  limits  of  the 
two  circles  generated  by  the  maximum  and  minimum  drift  errors.     See 
Figure  3. 

The  total  probable  error  of  position  can  now  be  determined.     This 
error,   c,   is  determined  by  the  equation: 

c  =  Vde2  +  X2  +  Y2 

13Ibid.,   p.   6-7. 
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Although  the  National  Search  and  Rescue  Manual  does  not  present  the 
development  of  this  equation  from  the  theoretical  background,  a  separate 
study      conducted  in  1964  derived  it  mathematically.     The  equation  was 
derived  using  the  Law  of  Cosines  and  the  root-mean-square  of  maxima 
and  minima  treating  the  three  error  vectors,  two  at  a  time. 


FIGURE  3 


GRAPHICAL  SOLUTION  OF  DRIFT  ERROR  (MINIMAX) 


U 


Furthermore,  if  the  values  for  X,  Y  and  de  are  considered  to  represent 
the  standard  deviation  for  a  normal  distribution  associated  with  each  of 
these  errors,  the  variance  would  be  X2,  Y2  and  de2  respectively.    Then  c 

Robert  L.  Armacost  and  Norman  T..  Saunders,    "Computer  Program 
of  USCG  Search  and  Rescue  Procedures"  (unpublished  Bachelor's  research 
paper,  The  United  States  Coast  Guard  Academy,  New  London,    1964). 

unnumbered. 


26 


would  represent  the  standard  deviation  of  a  normal  distribution  representing 
.  the  sum  of  these  errors.     The  1964  study  of  Armacost  and  Saunders  also 
reached  this  conclusion.  ^ 

The  American  Practical  Navigator^  devotes  a  chapter  to  the 
consideration  and  study  of  navigational  errors.    Among  the  errors  included 
in  this  discussion  are  constant  errors,   random  errors,,   periodic  errors,  and 
systematic  errors. 

Then  attention  is  directed  to  the  combination  of  these  errors  in  any 
given  line  of  position  or  navigational  fix.     The  distribution  of  these  various 
errors  and  the  application  of  probability  of  error  is  applied.     The  resulting 
conclusion  is,    "If  a  number  of  random  errors  are  combined,  the  result 
tends  to  follow  a  normal  curve  regardless  of  the  shape  of  the  individual 

errors "*■'    Furthermore,  applying  statistical  theory  the  authors  state, 

,  ".,.  .the  standard  deviation  of  all  the  errors  is  found  by  squaring  each 
individual  error,  adding  the  results,  and  taking  the  square  root  of  the  sum.  *?■" 
They  further  point  out  that  a  given  line  of  position,  as  finally  determined 

and  plotted,  "...  may  include  30  errors  or  more.  "^ 

.' 

15 Ibid. 

"Nathaniel  Bowditch  and  others,  American  Practical  Navigator. 
H.O.   Pub.  No.  9,    (1958  edition,  Washington:  U.  S.  Navy  Hydrographic 
Office,    1958),   pp.  678-688. 

17Ibid. ,   p.  683. 

18Ibid. 

19 Ibid. ,   p.   684. 
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Although  Dutton,   in  his  book,   Navigation  and  Nautical  Astronomy, 
does  not  comment  directly  upon  the  distribution  of  navigational  errors, 
he  does  discuss  celestial  fixes  for  surface  craft  and  the  magnitude  of  the 
errors  encountered*     He  states  that,    ".  , .  under  average  conditions  an 
experienced  navigator  can  obtain  a  fix  that  is  seldom  in  error  by  more 
than  a  mile.  "20    For  celestial  navigation  in  aircraft,    "Under  good 
conditions.  .  .a  fix  accurate  to  5  to  10  miles  is  considered  good. . .  "  but 
"...  individual  sights  taken  by  experienced  navigators  may  be  more  than 
100  miles  in  error  in  extreme  cases,  "^l 

In  considering  loran  fixes  he  states: 

"In  general,   ground  waves  produce  a  line  of  position 
accurate  to  1.  5  miles  or  better  over  80%  of  the  area 
covered  by  this  stations.     Sky  waves  produce  a  maximum 
error  of  5  to  7  miles  for  80%  of  the  coverage  area»     The 
accuracy  of  a  position  is  increased,  of  course,   if  three 
or  more  lines  of  position  are  used  and  the  various  lines 
weighted,  by  their  position  relative  to  the  stations,  type 
of  wave,  angle  of  cutting  other  lines,   etc.  "22 

Having  determined  the  total  probable  error  of  position,   the  next  step 

is  to  use  this  information  in  determining  the  limits  of  the  area  to  be 

searched.     Initially,  the  search  area  must  be  sufficiently  large  to  insure 

that  the  total  probable  error  of  position  is  included.     The  theoretically 

ideal  search  area  would  be  a  circle  centered  at  datum  and  with  a  search 

20Benjamin  Dutton,   Navigation  and  Nautical  Astronomy  (tenth 
edition,  Annapolis:  United  States  Naval  Institute,    1951),    p.   332. 
■■■■■ 
on  ■ 

2*Ibid. ,   p.   537, 

22 Ibid.,    p.   228o 


radius  equal  to  or  greater  than  the  total  probable  error  of  position,   c. 

The  search  area  is  expanded  for  each  subsequent  search  following 
an  unsuccessful  search  by  increasing  the  search  radius,   R.     The  number 
of  searches  conducted  is  not  limited  to  a  specific  number  since  adverse 
weather  conditions  at  the  scene  may  reduce  the  effectiveness  of  each 
sweep.     Such  conditions  would,   therefore,  necessitate  additional 
searches  over  and  above  the  number  used  under  good  search  conditions. 

Each  succeeding  search  will  be  centered  on  the  datum  point  at  the 
time  the  particular  search  is  to  be  started.    This  means  that  the  area  where 
the  survivors  are  most  probably  located,  namely  datum,  will  be  searched 
more  frequently  than  the  outer  limits  of  the  search  area,,     The  outer  limits 
are  extended  for  each  succeeding  search.     The  search  radius,  which 
determines  the  outer  limits,   is  found  by  multiplying  the  total  probable 
error  of  position,   c,   by  a  safety  factor.    The  safety  factor  is  increased 
for  each  search  to  reflect  the  increased  total  probable  error  of  position 
which  has  resulted  from  increased  drift  and  the  associated  increased  drift 
error.     The  values  used  in  determining  the  search  radius  are  as  follows: 23 

Search  Safety  Factor  Search  Radius 

1  1.1  1.  1  x  c 

2  1.  6  1.  6  x  c 

3  2.0  2.0  x  c 

4  2.3  2.3  x  c 

5  2.  5  2.  5  x  c 


National  Search  and  Rescue  Manual,   op.  cit. ,   p0   6-10< 
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The  time  required  to  complete  a  search  will  be  a  function  of  many 
things,   such  as  the  number  of  search  units  participating,   the  speeds 
used  by  the  search  units,  the  type  of  pattern  used  in  the  search,  and 
the  size  of  the  area  to  be  searched.     If  several  hours  elapse  in  conducting 
a  search,   it  may  be  necessary  to  completely  recompute  the  total  probable 
error  of  position,   qt  ,to  reflect  the  increased  drift  error. 

As  indicated  earlier,  the  theoretical  search  area  would  be  a  circle. 
However,,  to  attempt  to  search  while  navigating  within  the  limits  of  a 
circle  is  not  particularly  suitable  for  many  of  the  different  types  of 
search  patterns  commonly  used.     Therefore,   the  search  area  is  usually 
squared  off  to  provide  a  rectrangular  search  area  with  a  side  of  length 
2R  and  a  total  area  equal  to  4R^.     This  search  area  provides  a  more 
suitably  defined  area  for  navigational  purposes. 

Figures  4  and  5  show  the  concept  of  the  "repeated  expansion"  search 
concept  for  a  stationary  datum  and  a  moving  datum.     This  concept  calls 
for  five  searches  following  consecutively  using  the  search  radii  indicated 
previously.    This  provides  for  the  greatest  concentration  of  search  efforts 
in  the  area  where  the  survivors  are  most  probably  located  and  minimum 
efforts  in  the  areas  of  least  probable  location,,  ^4 


24Ibid. ,   p.   6-12 
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FIGURE  4- 
REPEATED  EXPANSION  SEARCH  AREAS  FOR  STATIONARY  DATUM  POINT 
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FIGURE  5 
REPEATED  EXPANSION  SEARCH  AREAS  FOR  MOVING  DATUM  POINT 
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Having  determined  the  datum  point  and  search  area,   the  RCC  Controller 
must  then  determine  what  search  units  are  available,   how  long  it  will  take 
them  to  reach  the  scene,  what  type  of  search  pattern  is  most  appropriate, 
and  how  much  time  is  available  for  searching  (limiting  factors  include 
aircraft  fuel  and  hours  of  daylight  remaining  before  dark).    Considering 
these  factors,  the  RCC  Controller  must  make  a  value  judgment  of  the 
search  area  calculations  and  perhaps  modify  the  area  by  altering  the  size 
„of  the  search  radius  or  the  track  spacing  of  the  search  units  or  botho     This 
is  known  as  a  "feasibility  check,  "25    This  is  a  decision  which  must  be 
made  on  the  basis  of  the  specific  situation  which  exists  for  each  individual 
distress  incident,, 

II.       STUDIES  DEALING  WITH  PARACHUTE  DRIFT  AND  OCEAN  CURRENT  DRIFT 

Parachute  drift. 


In  1964,  while  at  the  United  States  Coast  Guard  Academy,  Cadets, 
First  Class,  Armacost  and  Saunders,   conducted  a  study26  of  the  general 
problems  associated  with  search  and  rescue  operations.    A  major  portion 
of  their  efforts  was  directed  in  examining  the  problem  of  free-fall  and 
parachute  drift  of  survivors  who  were  required  to  bail  out  of  a  distressed 
aircraft, 

25Ibid^,    p0   7-11. 
'"Armacost  and  Saunders,   oj).  cit. 
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Their  examination  of  the  table  provided  in  the  National  Search  and 
Rescue  Manual  lead  them  to  believe  that  the  values  contained  therein  were 
related  linearly.     They  also  concluded  that  the  table  was  based  on  a  model 
where  the  survivor's  initial  velocity  at  time  of  bail-out  was  zero  in  both 
the  horizontal  and  vertical  planes.    Also,   the  table  utilizes  an  average 
wind  between  bail-out  altitude  and  the  surface.     It  makes  no  provision  for 
the  wide  variance  in  both  wind  velocity  and  direction  which  may  be 
encountered  at  different  altitudes. 

Accordingly  they  set  out  to  derive  equations  which  would  relate 
the  force  of  gravity,   effect  of  winds  in  a  horizontal  direction,   motion  due 
to  the  velocity  of  the  aircraft,,   and  air  resistance.     They  did  this  by 
developing  equations  for  the  vertical  component  of  free  fall,   the  horizontal 
component  of  motion  due  to  wind,   the  horizontal  component  of  motion  due 
to  an  initial  velocity  equ,al  to  that  of  the  speed  of  the  aircraft,  the  effect 
of  air  resistance  or  drag  as  a  function  of  air  density  which  in  turn  is  a 
function  of  barometric  pressure  and  temperature,   the  relationship  of 
kinematic  viscosity  of  air,   air  density  and  altitude,   and  finally,  the 
evaluation  of  the  various  constants  appearing  in  the  various  equations. 
Having  completed  this  extensive  development,  they  then  wrote  a  basic 
flow  diagram  for  their  program.     That  portion  of  their  work  dealing  with 
the  parachute  drift  problem  has  been  reproduced  in  its  entirety  in 
Appendix  C. 
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In  1965,  Cadet,   First  Class,   Darvy  M.  Cohan,  also  of  the  United 
States  Coast  Guard  Academy,   continued  the  work  started  by  Armacost  and 
Saunders.     Cohan  refined  and  developed  a  computer  program^'  incorporating 
the  results  of  Armacost  and  Saunders;  but,   in  addition,   he  investigated  the 
effects  of  changing  wind  direction  and  velocity  using  three  different 
altitudes:  (1)  surface  winds,    (2)  medium  altitude  winds,   and  (3)  high 
altitude  winds.     These  winds  were  selected  since  information  on  their 
magnitudes  is  usually  available  from  various  agencies  providing 
meteorological  data. 

In  addition,   he  considered  the  effects  of  three  possible  parachute 
opening  conditions  which  included  the  two  possible  extremes,   namely  (l) 
that  the  parachute  opened  immediately  upon  leaving  the  aircraft,  or  (2) 
the  parachute  did  not  open  at  all,  and  the  possibility  of  opening 
automatically  at  14,  000  feet.     This  latter  altitude  is  the  one  at  which  most 
military  parachutes  are  designed  to  open  using  a  barometric  automatic 
release  mechanism. 

Cohan  employed  the  Euler  method  of    integration  in  his  program. 
He  tested  his  program  on  an  IBM  1620  computer  and  found  that  the  running 
time  was  slightly  less  than  ten  minutes.    Accordingly,   he  recommended 

97 

Darvy  M.  Cohan,    "Free  Fall  and  Parachute  Drift"  (unpublished 
Bachelor's  research  paper,   The  United  States  Coast  Guard  Academy,  New 
London,    1965). 
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that  future  work  in  this  area  should  include  consideration  of  utilizing  more 
sophisticated  techniques  such  as  the  Runge-Kutta  method.     Furthermore,   he 
pointed  up  the  weakness  that  exists  in  that  no  consideration  has  been  given 
to  the  effect  of  turbulent  weather  conditions  with  the  associated  vertical 
shear,   updrafts,  and  downdrafts  which  may  have  contributed  to  precipitation 
of  the  distress  itself. 

Appendix  D  contains  the  basic  equations,   flow  diagrams,  and  computer 
program  for  the  IBM  1620  which  Cohan  developed. 

Ocean  current  drift. 

As  military  and  commercial  maritime  operations  continue  to  grow, 
the  need  for  more  and  more  accurate  information  dealing  with  surface  (and 
sub-surface)  currents --their  prediction  and  behavior --increases  in 
importance.     Oceanographers  and  meteorologists  have  long  been  studying 
this  problem*     Many  theories  have  been  advanced;  however,  obtaining 
sufficient  experimental  data  to  test  these  theories  has  proved  to  be  a 
major  obstacle. 

The  Fleet  Numerical  Weather  Facility  (FNWF)  at  Monterey, 
California,  has  recently  developed  a  computerized  system  for  producing 
synoptic  analyses  of  surface  currents  based  on  the  distribution  of  sea 
surface  temperatures  and  the  climatological  temperature  field  at  a  depth 
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of  200  meters.    Commander  Wu   Ec  Hubert,   USN,   presented  a  paper28  at 
the  Institute  of  Navigation  in  San  Francisco  in  December  1964  is  which 
he  discussed  the  general  method  used  in  the  Fleet  Numerical  Weather 
Facility's  program.    This  program  is  designed  to  evaluate  the  transport 
due  to  the  permanent  flow  component  and  the  component  due  to  wind  and 
waves  thus  producing  the  overall  pattern  of  surface  currents.    Commander 
Hubert's  paper  appears  in  Appendix  E. 

Taivo  Laevastu  has  conducted  an  extensive  review  of  the  materials 
dealing  with  surface  currents  in  connection  with  his  work  at  the  Fleet 
Numerical  Weather  Facility  and  his  preparation  of  the  book  Qceanoqraphic 
Forecasting  for  the  Navy  Oceanographic  Office,     In  Chapter  Ten,   Doctor 
Laevastu  examines  in  detail  the  various  aspects  of  "Prediction  of  Surface 
Currents.  "^   Among  the  factors  investigated,  he  discusses  the  evaluation 
of  permanent  flow,  wind  currents,,   mass  transport  by  waves,  inertia  currents, 
tidal  and  "hydraulic  currents,  "  currents  caused  by  the  change  in  atmospheric 
pressure,  and  other  current  components.  ^0    He  further  discusses  the 


28 
W.  E.  Hubert,   "Computer  Produced  Synoptic  Analyses  of  Surface 

Currents  and  Their  Application  for  Navigation"  (paper  read  at  the  1964 

National  Marine  Navigation  Meeting,  Institute  of  Navigation,  San 

Francisco,   December  7-8,    1964). 
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Taivo  Laevastu,   "Prediction  of  Surface  Currents,  "  Oceanographic 

Forecasting  (a  book  presently  under  preparation  for  U.  S.  Naval 

Oceanographic  Office  under  Contract  Number  N62306-1129). 


30Ibid_.,   pp.    10-1  -  10=26. 
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.evaluation  of  change  of  current  speed  and  direction  with  depth,  the 
influence  of  water  depth,  configuration  of  coast  and  islands,  nearshore 
circulation,   current  boundaries  and  eddies,  and  deep  currents.  31 

Further  studies  relating  the  work  of  the  Fleet  Numerical  Weather 
Facility  and  the  problems  associated  with  surface  weather  phenomena 
are  found  in  their  Technical  Notes  Numbers  Five ^ 2  and  Eighths  published 
•in  February  1965.    These  Technical  Notes  discuss  the  sea-air  interactions 
on  a  synoptic  scale  and  the  analyses  and  forecasting  of  sea  surface 
temperatures,  both  of  which  are  directly  related  to  the  paper  presented 
by  Commander  Hubert  in  December  1964  and  referred  to  earlier. 

The  U.  S.  Naval  Oceanographic  Office  has  been  vitally  interested 
in  the  problems  associated  with  ocean  current  prediction.   Although 
military  operations  involving  anti-submarine  warfare  have  given  impetus 
to  studies  in  this  area,  the  problems  associated  with  prediction  of  the 
location  of  freely  drifting  mines  and  rafts  and  the  associated  areas  for 
searching  for  these  objects  have  also  been  under  investigation. 

Donald  A.   Burns  of  the  U.  S.  Naval  Oceanographic  Office  has 


31  Ibid. ,   pp.   10-23  -  10-32. 

32W.  E.  Hubert,   "U.  S.  Fleet  Numerical  Weather  Facility  Activities 
Relating  to  Sea-Air  Interactions  on  a  Synoptic  Scale"  (Technical  Note 
Number  Five.     Monterey:  Fleet  Numerical  Weather  Facility,    1965). 

°P.   M.  Wolff,  L.  P.  Carstensen  and  T.  Laevastu,   "Analyses  and 
Forecasting  of  Sea  Surface  Temperature  (SST)"  (Technical  Note  Number 
Eight.     Monterey:  Fleet  Numerical  Weather  Facility,    1965). 
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developed  a  Fortran  Computer  Program  which  defines  a  drift  ellipse  pattern 
for  the  free-floating  object.  ^4  The  essentials  of  this  drift  ellipse  program 
are  described  in  A  Statistical  Rose  Program  developed  by  Walter  E„  Yergen 
and  published  by  the  Oceanographic  Office.  ^^ 

Basically,   Bums'  program  analyzes  input  data  consisting  of 
"grouped  frequency  current  observations.  "36    This  analysis  produces 
various  parameters  which  are  then  plotted  to  outline  that  area  in  which 
recovery  is  most  probable.    This  is  accomplished  by  determining  vector 
means,  vector  variances  and  covariances,  determining  the  inclination  of 
the  principle  axis  of  the  distribution,  and  the  standard  deviations  along 
these  principle  axis.    Next,  within  specified  limits  for  speed  and  direction 
classes,  the  frequencies  for  the  normal  surface  are  integrated.  37    This 
Drift  Ellipse  Program  is  found  in  Appendix  F. 

III.      NOTE  ON  ADDITIONAL  LITERATURE  RELATING  TO  THE 
OVERALL  SEARCH  AND  RESCUE  PROBLEM 

In  undertaking  this  project  it  was  necessary  that  a  thorough  general 

34Donald  A.   Burns,    "Instructions  for  Surface  Drift  Ellipse  Program" 
(Washington:  U.  S.  Naval  Oceanographic  Office,  June  1962). 

35 

Walter  E.  Yergen,  A  Statistical  Rose  Program,  SP-64  (Washington: 

Evaluation  Branch,  Oceanographic  Analysis  Division,   Marine  Sciences 

Department,  U.  S.  Naval  Oceanographic  Office,  October  1962). 

Burns,  op.  cit. ,   passim. 

3 'Yergen,  op.  cit. ,   passim. 
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understanding  of  all  aspects  of  the  search  and  rescue  problem  be  acquired 
before  an  attempt  was  made  to  sub-optimize  that  portion  which  was 
selected.    Accordingly,  the  Bibliography  at  the  end  of  this  paper  is 
annotated  to  reflect  all  the  pertinent  literature  which  was  reviewed  in 
the  process  of  this  study. 
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CHAPTER  III 
THE  STUDY 

I.       METHOD 

The  method  of  research  employed  in  this  study  was  basically  that  of 
the  development  of  a  conceptual  model  designed  for  a  computer  solution 
expressing  the  doctrine  outlined  in  the  National  Search  and  Rescue  Manual 
for  the  determination  of  a  search  area.     Through  the  application  of 
mathematical  techniques  relationships  of  the  numerous  variables  were 
developed.    These  relationships  were  then  outlined  in  a  series  of  computer 
flow  diagrams,  for  various  subroutines,  and  the  computer  programs  were 
written  in  Fortran  60  computer  language.    After  each  subroutine  was 
thoroughly  checked  and  tested,  the  final  program  was  assembled,  tested 
and  evaluated  using  hypothetical  distress  cases. 

II.      PROCEDURE 

The  problem  of  determining  a  search  area  for  any  given  distress 
incident  involves  three  basic  steps: 

1.  Determine  the  datum  point  at  any  specified  time. 

2.  Determine  the  probable  position  errors  associated  with  the  datum 
point. 

3.  Determine  the  limits  of  the  area  which  will  be  included  in  the 
search. 
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Each  of  these  steps  was  examined  individually  and  in  detail.    The 
examinations  were  accomplished  incrementally  by  determining  what 
output  was  required,  what  were  the  inputs  available  or  necessary  to 
develop  the  desired  output,  and  how  were  the  inputs  and  output  related. 
The  limitations  outlined  in  Chapter  I  were  carefully  observed. 


START 


DETERMINE 
DATUM 


<- 


i 


DETERMINE 
PROBABLE 
POSITION  ERROR 


<J- 


SL 


DETERMINE 
SEARCH 
AREA 


VARIOUS, 
INPUTS 


VARIOUS/ 

INPUTS> 


VARIOUS; 
INPUTS 


xf 


FIGURE  6 
THE  BASIC  PROCESS  OF  DETERMINING  THE  SEARCH  AREA 

The  following  sectionspf  this  chapter  will  treat  each  of  the  above  areas 
in  order.     Program  flow  diagrams  will  reflect  the  general  logical  sequence 
which  was  employed  although  the  terminology  used  will  be  broad  enough 
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to  permit  the  development  of  similar  programs  using  computer  languages 
other  than  Fortran  60.     Standard  symbols  will  be  used  throughout  the 
diagrams„    An  index  to  these  symbols  will  be  found  in  Appendix  A.     The 
detailed  computer  program  in  Fortran  60,   complete  with  definitions  of  the 
applicable  variable  names  used,  will  be  found  in  Appendix  B. 

Determination  of  datum  at  any  specified  time,  _tj[. 

When  a  distress  case  has  entered  into  the  "Distress"  phase,  the 
RCC  Controller's  first  task  is  to  determine  the  datum  point.     Since  this 
datum  point  will  be  continuously  moving  as  time  passes,   he  must  select 
the  specified  time,  ti,   for  which  datum  is  to  be  determined.     The  initial 
reported  position  of  the  distress  incident  will  be  datum  t0  where  t0  is 
the  time  at  which  the  distress  incident  occurred.     The  time  selected,  ti, 
may  be  the  estimated  time  of  arrival  of  search  units  at  the  scene,  or  it 
may  be  the  predicted  position  at  any  other  time  during  the  period  of  the 
case  duration.     It  will  be  redetermined  during  each  search  and  used  as 
the  center  of  subsequent  searches.    What  then  is  the  information  needed 
by  the  RCC  Controller  to  determine  datum  ti? 

The  basic  information  needed  will  be  the  initial  distress  incident 
position,   datum  t0,  or  the  previous  datum  if  this  is  a  redetermination 
of  datum,  the  time  of  the  previous  datum,  the  time  for  which  the  new 
datum  is  being  determined,  and  the  drift  vectors  which  have  been  acting 
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upon  the  previous  datum  during  the  time  interval  which  will  have  elapsed. 
See  Figure  7 -A. 
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THE  FACTORS  ASSOCIATED  WITH  THE  FORECASTED  DATUM 

The  following  specific  information  is  necessary: 

1.  The  type  of  unit  in  distress,  i.e.,    an  aircraft,  a  ship,  or  a  small 
craft  (this  latter  portion  was  not  developed  in  this  study). 

2.  The  initial  distress  incident  position.  This  position  will  be  expressed 
in  terms  of  latitude  and  longitude. 
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3.  The  time  associated  with  the  initial  distress  incident  position  in 
terms  of  the  date,  hour  and  minute  expressed  in  Greenwich  Mean  Time 
(GMT). 

4.  Based  on  the  type  of  unit  in  distress  and  a  knowledge  of  the  type 
of  survival  equipment  carried  on  that  type  of  distressed  unit,  the  RCC 
Controller  must  specify  what  type  of  survival  craft  is  believed  to  be  the 
object  of  the  search.     This  information  may  have  been  furnished  in  the 
message  reporting  the  distress  incident. 

5.  Weather  information  at  the  scene  of  the  distress  incident  including: 

a.  The  present  surface  wind  velocity  expressed  in  knots  and  the 
direction  from  which  it  is  blowing  in  degrees  from  true  North. 

b.  The  surface  winds  for  the  twenty-four  hour  period  preceding 
the  time  for  which  the  datum  is  being  determined.    The  velocity  and 
directions  are  needed  and  must  be  expressed  in  knots  and  degrees  from 
true  North. 

.    6.  The  average  sea  current  set  and  drift  expressed  in  knots  and 
degrees  from  true  North.    This  information  is  extracted  from  the  Atlas  of 
Surface  Currents  for  the  appropriate  area  and  month.    At  the  same  time 
the  steadiness  of  this  current  is  also  determined  and  noted. 

7.  The  specific  time,  tj,   for  which  datum  is  to  be  determined.    This 
will  be  expressed  in  terms  of  4a te,   hour,  and  minute  using  GMT. 

Using  these  specific  inputs,  the  drift  components  of  leeway,  local 
wind  current,  and  average  sea  current  are  determined  as  follows: 
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Leeway.    Since  leeway  is  the  motion  of  the  survival  craft  resulting 
from  the  action  of  the  surface  wind  upon  the  craft,  the  first  factor  to  be 
considered  is  the  type  of  survival  craft  believed  to  be  employed.     If  the 
survivors  are  believed  to  be  floating  in  the  water,   either  unaided  or  with 
life  jackets,   the  effect  of  the  local  surface  winds  is  assumed  negligible. 
In  this  case  the  leeway  component  becomes  zero  and  the  computations 
continue  to  the  next  part  of  the  program. 

If,   however,  the  survival  craft  is  believed  to  be  either  a  raft  or  a 
boat,   a  leeway  component  will  exist  and  must  be  determined.     Considering 
the  case  of  a  raft  first,   its  motion  will  be  directly  downwind.     Therefore, 
the  input  of  the  local  surface  wind  direction  plus  180  degrees  will  yield 
the  direction  of  the  leeway  component.     The  rate  at  which  the  raft  will 
drift  is  a  function  of  the  local  surface  wind  velocity  and  also  depends 
upon  whether  or  not  the  survivors  in  the  raft  stream  a  drogue  to  retard 
their  motion  in  leaving  the  site  of  the  distress  incident.     Therefore,  two 
possible  drift  rates  must  be  considered;  (1)  the  rate  of  drift  with  a  drogue, 
and  (2)  the  rate  of  drift  without  a  drogue.     This  results  in  the  determination  of 
a  maximum  value  and  a  minimum  value  of  the  leeway  component.    When  the 
leeway  components  (in  knots)  are  multiplied  by  the  elapsed  time,  At,   the 
total  maximum  and  minimum  leeway  components  are  determined. 

The  treatment  in  the  case  of  a  boat  is  similar  to  that  of  a  raft. 
However,  the  direction  of  drift  of  a  boat  has  been  found  to  vary  as  much 
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as  40°  on  either  side  of  the  downwind  direction.     Therefore,   the  maximum 

and  minimum  leeway  vectors  must  also  include  maximum  and  minimum 

vectors  40°  on  either  side  of  the  downwind  direction  to  reflect  this 

divergence.     This  provides  a  total  of  six  leeway  vectors  of  which  the 

two  downwind  vectors  are  used  in  computing  datum  while  all  six  vectors 

are  considered  in  the  determination  of  the  drift  error,, 

In  this  specific  study,  the  magnitude  of  the  minimum  leeway  component 

was  determined  by  one  of  two  equations  based  on  leeway  curves  given  in 

the  National  Search  and  Rescue  Manual.  38   Analysis  of  the  curves  showed 

that  leeway  for  a  raft  with  a  drogue  could  be  closely  approximated  by  two 

straight  lines.     The  first  line  approximated  leeway  for  wind  velocities 

equal  to  or  less  than  five  knots: 

Leeway  in    n.  miles    =    (wind  velocity  in  knots  x  0.  8) 
hr.  24 

The  second  line  approximated  leeway  for  surface  wind  velocities  greater 

than  five  knots; 

Leeway  in    n.  miles    =    (wind  velocity  in  knots  +2.5) 
hr.  1.87  x  24 

The  magnitude  of  the  maximum  leeway  component  (no  drogue  streamed) 

was  also  determined  from  the  curves.     However,  the  leeway  curve  values 

had  been  stored  in  memory  and  were  furnished  directly  as  a  function  of  the 

wind  velocity.     Figue  7-B  is  the  flow  chart  for  the  program  used  in 

determining  the  leeway  components. 


38 

National  Search  and  Rescue  Manual,  op.  cit. ,    p.   6-4. 
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FIGURE  7-B 
FLOW  DIAGRAM  FOR  DETERMINING  LEEWAY  DRIFT  COMPONENTS 
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Applying  the  elapsed  time  to  the  leeway  magnitudes  and  then 
combining  this  product  with  leeway  direction  yields  the  leeway  vectors. 

Local  wind  current «,     The  local  wind  current  is  treated  in  the  manner 
^described  in  Chapter  II.     However,   the  mean  wind  direction  and  velocity 
are  determined  manually  with  the  number  of  calculations  involved  depending 
upon  the  extent  of  the  surface  wind  variations  during  the  preceding  twenty- 
four  hour  period.     The  mean  wind  for  the  preceding  twenty-four  period  is 
then  used  as  an  input  for  the  computer. 

The  magnitude  of  the  current  generated  by  the  local  wind  is  a  function 

of  the  mean  wind  velocity.     Examination  of  the  table  relating  the  velocity 

of  the  wind  current  to  the  mean  wind  speed  in  the  National  Search  and 

Rescue  Manual  suggested  a  linear  relationship.         The  data  from  the  table 

were  plotted  and  the  following  equation  provided  the  relationship: 

Velocity  of  the  local  wind  current    =    (mean  wind/1.  18) 
(in  knots)  24 

The  latitude  of  the  initial  distress  position  (or  previous  datum)  was 

compared  to  the  10°N  latitude  criterion  to  determine  the  divergence  of  the 

local  wind  current  resulting  from  the  (3oriolis  force.     If  the  distress  inptdent 

latitude  was  less  than  or  equal  to  10°N,  the  direction  of  the  local  wind 

current  became  the  direction  of  the  mean  local  wind  +  180°.     If  the  distress 

incident  latitude  was  greater  than  10°N,  the  direction  of  the  local  wind 

39Ibid. 
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current  diverged  30°  to  the  right  of  the  downwind  direction  of  the  mean 
local  wind. 

By  applying  the  elapsed  time  to  the  local  wind  current  velocity  and 
then  combining  this  with  the  direction  of  the  local  wind  current/  the 
final  local  wind  current  drift  vector  was  determined.     See  Figure  7-C. 

Average  sea  current.     This  component  is  simply  the  values  determined 
from  the  Atlas  of  Surface  Currents  for  set  and  drift  with  the  latter  portion 
multiplied  by  the  elapsed  time.     If  the  degree  of  steadiness  was  less  than 
35%,   this  component  is  assigned  a  zero  value  by  the  RCC  Controller. 
See  Figure  7-D. 

Resultant  total  drift.      The  final  resultant  total  drift  vector  for  the 
elasped  time  period  is  the  vector  sum  of  the  leeway  components,  the 
local  wind  current  component,  and  the  average  sea  current  component., 
Since  this  involves  maximum  and  minimum  components,   maximum  and 
minimum  total  drift  vectors  will  be  generated.     In  accordance  with  the 
assumptions  stated  in  Chapter  I,   the  new  datum  will  be  found  at  a  point 
midway  between  the  maximum  drift  position  and  the  minimum  drift  position 
found  by  using  the  downwind  leeway  vectors.     See  Figure  2  in  Chapter  I, 
page  11,  for  the  general  drift  picture.     Figure  8  shows  the  drift  picture 
for  a  life  raft,  and  Figure  9  shows  the  drift  picture  for  a  boat  with  the 
associated  divergence  of  the  leeway  vectors.     In  the  latter  case  the  most 
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probable  position  of  the  survivors  is  assumed  to  lie  midway  between  the 
maximum  and  minimum  positions  on  the  downwind  vector.     However,  the 
extreme  drift  positions  associated  with  the  divergence  of  the  leeway 
vectors  serve  as  input  data  for  a  portion  of  the  problem  phase  relating 
to  position  errors. 

In  the  specific  program  developed  by  this  study  the  problem  of 
vector  addition  was  treated  in  two  steps.     The  first  step  was  to  convert  the 
directions  of  the  vectors,  which  were  expressed  in  degrees  measured 
clockwise  from  True  North,   to  the  standard  Cartesian  coordinate  system 
used  by  the  computer.     Computer  angles  are  measured  in  radians.     The 
conversion  was  accomplished  by  the  relationship: 

%  =  *f-  -  Or 

Where  #r  is  the  angle  expressed  in  radians  as  used  in  the  computer 
calculations,  and  0^  is  the  direction  of  the  component  converted  from 
degrees  True  into  radians. 

The  second  step  in  treating  the  vector  addition  problem  was  to  convert 
the  various  components  of  leeway,   local  wind  current,  and  average  sea 
current  into  their  respective  x  and  y  components.     The  final  resultant 
vector  is  obtained  by  summing  the  x  components,   summing  the  y  components, 
and  then  applying  the  Pythagorean  Theorem  to  find  the  resultant  vector. 

Actually,   the  various  drift  vectors  were  broken  into  the  x  and  y 
components  and  carried  as  subscripted  variables  throughout  the  program 
and  applied  where  needed. 
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Determination  of  total  probable  error  of  position. 

The  total  probable  error  of  position  is  a  function  of  three  basic  errors: 
(1)  the  initial  error  associated  with  the  distress  incident  position   (2)  the 
navigation  error  associated  with  the  search  units ,  and  (3)  the  error 
associated  with  the  determination  of  the  drift  components.    See  Figure 
10-A. 
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FIGURE  10-A 
COMPONENTS  OF  TOTAL  PROBABLE  ERROR  OF  POSITION 

The  following  specific  information  is  needed  to  determine  the  values 
of  these  errors: 

1.  The  type  of  unit  in  distress. 

2.  The  type  of  search  unit  which  will  be  dispatched. 

3.  The  classification  of  the  method  of  navigation  employed  to  determine 
the  initial  distress  incident  position.    If  the  method  of  navigation  is 

unclassified  or  is  classified  as  dead  reckoning  (DR),  the  last  known  fix 
of  the  distress  unit  must  be  ascertained. 
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4.  The  classification  of  the  method  of  navigation  which  will  be 
employed  by  the  search  unit  in  proceeding  to  the  scene  of  the  distress 
incidento 

5.  The  present  position  of  the  search  unit  must  be  determined.     This 
present  position  of  the  search  unit  is  treated  as  a  fix.     It  must  be  expressed 
in  terms  of  latitude  and  longitude. 

6.  The  magnitude  of  the  maximum  drift  vector0 

7.  The  magnitude  of  the  minimum  drift  vector. 

Applying  these  specific  inputs,   the  three  types  of  errors  are  determined 
as  follows: 

Initial  position  error  of  the  distress  incident.     This  error  is  strictly 
a  function  of  the  type  of  craft  in  distress  and  the  method  of  navigation  used 
to  determine  the  initial  position  of  the  distress  incident.     In  a  general 
program  the  procedure  is  to  test  for  the  type  of  craft  in  distress  which  then 
determines  which  portion  of  the  program  is  to  be  used.     The  next  step  is 
to  determine  the  classification  of  the  method  used  in  determining  the 
reported  incident  position.     Using  the  values  for  initial  position  error 
cited  in  Chapter  II  for  the  appropriate  navigational  method  produces  the 
desired  initial  position  error,  X. 

If  the  classification  of  the  navigational  method  is  that  of  a  position 
based  upon  radio  direction-finder  bearings  (DF),   the  type  of  DF  system 
used  must  be  further  identified.     The  two  types  of  DF  systems  commonly 

55 


used  are  the  commercial -military  DF  service  nets  and  the  Federal 
Communications  Commission  (FCC)  nets.     Each  of  these  nets  will  furnish 
a  classification  of  their  position  fix  which  indicates  the  error  associated  with 
the  position.     These  classifications,   which  are  used  in  the  program,  are 
given  below.  40 

Fix  Classification        Associated  Position  Error  As  A  Radius  In  Nautical  Miles 

DF  Service  Net  FCCrDF  Net 

Class  A  Within  5  20  or  less 

Class  B  Within  20  40  or  less 

Class  C  Within  50  60  or  less 

Class  D  -  -  -  -  more  than  60 

Since  the  DF  classification  system  is  independent  of  the  type  of 
craft  involved,   one  portion  of  the  computer  program  handling  this  type  of 
position  classification  serves  for  all  other  branches  of  the  system 
regardless  of  the  type  of  craft  involved. 

Originally  the  specific  program  developed  by  this  study  followed  the 
general  program  shown  in  Figure  1Q-B.     The  various  navigational 
classifications  were  assigned  different  code  numbers  or  combinations  of 
codes.     These  codes  were  then  punched  into  the  input  data  cards.     Since 


40 Ibid. ,   p.   5-7  „ 
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it  was  obvious  that  the  operator  in  the  RCC  would  have  to  put  in  these 
classification  code  numbers  in  one  form  or  another,   it  was  decided  to 
change  the  code  numbers  to  the  actual  values  of  the  error  for  the  methods 
which  were  classified  as  a  Fix  or  DF  Fix„     This  limited  the  program 
calculations  to  those  involved  with  DR  classifications.     This  will  probably 
be  the  classification  most  generally  encountered  in  actual  practice. 

In  the  case  of  a  DR  classification,   the  error  is  equal  to  the  error 
associated  with  the  last  known  fix  plus  a  certain  percentage  of  the  distance 
traveled  from  that  last  fix  to  the  reported  DR  position.     The  percentage  of 
the  distance  used  is  dependent  upon  the  type  of  craft  involved.     The 
problem  encountered  in  this  determination  is  ascertaining  the  distance 
traveled  between  those  two  points. 

There  are  various  methods  available  for  determining  distances  between 
points  on  the  Earth's  surface.     The  method  selected  for  this  program  was 
that  of  Mercator  Sailing.     The  only  inputs  needed  for  the  computer  were 
the  latitudes  and  longitudes  of  the  two  points „     The  basic  equations  used 
in  solving  this  problem  were  those,  developed  by  Bowditch  in  his  book, 
American  Practical  Navigator.  41    They  are: 

D    -    c4c    .      and     tanC    =   22*.     , 


41 

Bowditch,   op.  cit. .   p.   227, 
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where:  D  =  the  distance  in  nautical  miles  between  1  and  2. 

X  -  the  difference  in  latitudes,   L-^  and  L~,   in  minutes. 

DLo  =  the  difference  in  longitudes,  X  -.  and  K    ,   in  minutes. 

m  =  the  difference  in  meridional  parts,    M-^  and  JVL,    (where 
meridional  parts  are  a  measure  of  the  arc  length  of  a 
meridian  between  a  given  latitude,   I4,   and  the  equator 
measured  on  a  Mercator  chart,   expressed  in  units  of 
1'  of  longitude  at  the  equator).  ^2 

M^  =  the  number  of  meridional  parts  between  the  equator  and 
latitude,   L^. 

The  values  of  M^  were  obtained  by  using  the  first  two  terms  of  the 

equation:  43 

Mi  =  7915,  704468  log  tan(45°  +  Lj/2)    -    23.  2689323^ 

-    0o  0525003^31^    -    0.000213sin5Li    -    ... 

Navigational  error  of  the  search  unit.    This  error  follows  the  same 
identical  pattern  of  solution  as  that  utilized  in  the  determination  of  the, 
initial  position  error,  X.     The  general  flow  diagram  is  shown  in  Figure 
10-C.    Only  the  position  of  the  last  known  fix  of  the  unit  concerned  will 
vary. 


42Ibid_. ,   p.  936, 

43ibid. ,   p„    1187.     "The  constants  used. .  .are  based  upon  Clarke's 
spheroid  of  1866.  .  .the  standard  reference  spheroid  used  for  charting 
North  America.  "    In  the  computer  program  only  the  first  six  significant 
digits  were  used  for  each  constant. 
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FIGURE   10-C 

FLOW  DIAGRAM  FOR  THE  DETERMINATION   OF   THE  NAVIGATIONAL 
ERROR  OF  THE   SEARCH   UNIT,    Y. 
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Drift  error.     The  methods  available  for  the  determination  of  drift 
error,   de,  were  outlined  in  Chapter  II.     The  computer  solution  for  de 
was  achieved  by  using  the  vector  summation  of  the  various  components 
of  drift  as  discussed  earlier.     These  vector  sums  produced  maximum  and 
minimum  drift  vectors.     In  the  case  of  a  raft  there  are  two  vectors,  one 
maximum  and  one  minimum.     For  a  boat  there  are  six  vectors,   three 
maximum  values  and  three  minimums0    These  vectors  are  tested  to  find 
the  positions  which  are  in  fact  farthest  from  and  closest  to  the  datum 
which  served  as  the  starting  point  for  the  drift  determination.     Using 
the  maximum  and  minimum  values  obtained  in  the    test  as  entering 
arguments,  a  table  44  stored  in  memory  produced  the  final  value  of  drift 
error,   de.     See  Figure  10-D. 

The  final  determination  of  the  total  probable  error  of  position,   c, 
is  then  accomplished  by  combining  the  three  separate  errors  in  accordance 
with  the  equation: 


c    =  Yx2    +    Y2    +    de2 
See  Figure  10-E. 

Determination  of  the  limits  of  the  search  area. 

Once  a  datum  point  has  been  established  for  a  specified  time,  ti, 


National  Search  and  Rescue  Manual,   op.   cit. ,    p.   6-9. 
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and  the  total  probable  error  of  the  position  has  been  determined,  the 
boundary  limits  of  the  search  area  may  be  defined.    These  boundaries 
may  be  expressed  in  either  of  two  ways.     The  first  is  in  the  form  of  the 
length  of  the  search  radius  which  extends  from  the  datum  point.    This 
method  is  used  for  vector  type  searches.    The  second  is  in  the  form  of 
geographical  coordinates  of  latitude  and  longitude  defining  a  rectangular 
search  area  centered  on  datum  with  a  length  and  width  each  equal  to 
twice  the  search  radius»    See  Figure  11» 
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FLOW  DIAGRAM  FOR  THE  DETERMINATION  OF  THE  SEARCH  AREA 
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The  specific  input  data  needed  to  define  the  search  area  is: 

1.  The  number  of  the  search  being  conducted.     This  is  needed  to 
ascertain  what  safety  factor  will  be  applied  in  the  calculation  of  the 
radius,,   R. 

2.  The  various  safety  factors  associated  with  each  different  search 
number.     These  are  the  factors  listed  in  Chapter  II,   page  29. 

3.  The  total  probable  error  of  position,   c.     This  has  been  computed 
in  the  preceding  phase. 

4.  The  coordinates  of  the  forecasted  datum  which  were  determined 
in  the  first  phase  of  the  overall  program. 

The  actual  computations  are  fairly  simple.     The  search  radius,   R, 
is  the  product  of  the  safety  factor  and  the  total  probable  error  of  position, 
c.     This  may  be  printed  out  along  with  the  coordinates  of  the  forecasted 
datum  to  provide  the  first  method  of  defining  the  search  area. 

In  the  second  method,  the  latitude  increment  will  be  equal  to  the 
search  radius  (in  nautical  miles)  divided  by  60  nautical  miles  per  degree 
of  latitude.     This  increment  of  latitude  is  added  to  the  forecasted  datum 
latitude  to  yield  the  northern  limit  and  subtracted  from  the  datum  latitude 
to  provide  the  southern  limit. 

In  determining  the  increment  of  longitude  equal  to  the  search  radius, 
the  datum  latitude  was  considered  to  be  the  middle  latitude,   Lm,  of  a 
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small  area  Mercator  plotting  sheet.    The  equation  relating  the  search 
radius  distance  in  miles  to  degrees  of  longitude  is^5 

Degrees  of  Longitude    =  ^ — -* R  miles 

*  60cosLm  miles/0  of  longitude 

This  increment  was  then  subtracted  from  the  datum  longitmde  to  obtain 
the  longitude  of  the  eastern  limit  and  was  added  to  the  datum  longitude 
for  the  western  limit. 

This  then  completes  the  basic  program  for  the  determination  of  the 
search  area  problem.     However,  the  specific  program  which  was  developed 
contains  three-  additional  items  of  interest.    The  first  two  of  these  are 
designed  to  provide  general  information  concerning  the  possible  conditions 
of  the  survivors.    This  information  assists  the  RCC  Controller  in  making 
the  difficult  decision  of  when  to  terminate  an  unsuccessful  search  and 
rescue  case.    It  will  also  provide  information  to  the  search  units 
-indicating  tentative  survival  times. 

The  first  of  these  two  subroutines  interrogates  a  table  stored  in  memory 
which  provides  information  expressing  the  effect  of  wind  and  air  temperatures 
©n  exposed  survivors.    This  table  converts  the  values  of  wind  and  air 
temperatures  at  the  scene  of  the  incident  to  the  equivalent  temperature 
on  exposed  flesh  if  the  wind  velocity  were  zero.  ^ 

45 

Bowditch,  op.  cit. ,   p.  89.    This  section  describes  the  various 

mathematical  relationships  between  latitudes,  longitudes,  and  distances 

on  small  area  Mercator  plotting  sheets. 

46 

This  information  was  obtained  from  a  "Wind  Chill  Chart"  developed 

by  the  U.S.  Army.    It  appeared  in  the  December  18,   1964  Commandant's 

Bulletin,  No.  51-64,  Supplement  No.  4.  published  by  the  U.  S.  Coast  Guard. 

This  table  was  reprinted  from  U.S.N.  Medical  Newsletter,   Vol.  32,  No.   12. 
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The  inputs  needed  are  two  in  number.    First,  the  air  temperature  at 
the  scene  is  needed  and  may  be  obtained  from  weather  reports  or  forecasts 
for  the  area.     The  second  input  is  that  of  the  velocity  of  the  local  surface 
winds  in  the  area.    This  has  been  entered  in  an  earlier  phase. 

These  two  inputs  serve  as  the  entering  arguments  for  the  table  which 
then  reads  out  the  desired  information  directly.    See  Figure  12. 
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FIGURE  12 

GENERAL  FLOW  DIAGRAM  FOR  EQUIVALENT  TEMPERATURES 
AT  VARIOUS  WIND  VELOCITIES 


The  second  subroutine  associated  with  the  condition  of  the  survivors 
deals  with  immersion  hypothermia,   that  is,  the  life  expectancy  of 
survivors  immersed  in  water  at  different  temperatures.     The  desired 
information  is  an  estimate  of  the  time  survivors  will  remain  conscious 
and  an  estimate  of  the  expected  time  they  have  before  death  becomes  an 
almost  certainty. 
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The  only  input  needed  for  this  portion  of  the  program  is  that  of  the 
sea  water  temperature  in  the  area  of  the  distress  incident.     This  information 
is  also  available  in  the  local  weather  reports  or  forecasts  for  the  distress 
area. 

This  input  information  is  the  entering  argument  into  tables  which 
provide  the  time  ranges  at  different  temperatures  for  consciousness  and 
time  until  death.  47    See  Figure  13, 

The  third  and  final  item  pertains  to  the  format  of  the  final  output. 

In  order  to  expedite  the  entire  process  of  transmitting  information  to  the 

search  units  which  may  be  assigned,  the  format  of  the  computer  printed 

output  was  designed  to  be  that  of  the  standard  Outgoing  Message  format. 

Instead  of  taking  the  information  from  a  given  computer  run,  transcribing 

this  information  to  message  blanks,  and  then  routing  it  to  the  communications 

center,   the  actual  print  out  serves  as  the  drafted  message.     The  only 

information  which  must  be  inserted  is  the  message  date-time  group  (DTG), 

address  or  name  of  search  unit  involved,  identification  of  the  distress 

incident,  and  the  releasing  officer's  approval,     This  also  eliminates  the 

possibility  of  transcription  errors  in  the  many  numbers  involved. 

Appendix  B.   contains  samples  of  the  print  out  format  for  the  first  and 

subsequent  searches  with  computer  calculated  figures  indicated  by 

i  I. 

47Climatological  and  Oceanograpnic  Atlas  for  Mariners  (Washington: 
Office  of  Climatology  and  Oceanographic  Analysis  Division,   Department  of 
Commerce,    1961),  Chart  No.    157. 

68 


underscoring  and  RCC  Controller  entered  information  indicated  by 
underscored  X's  (XXX).     In  the  space  after  the  last  standard  print  out 
paragraph  the  RCC  Controller  can  enter  any  additional  information  that 
he  deems  appropriate.  j 
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FLOW  DIAGRAM  FOR  IMMERSION  HYPOTHERMIA 
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CHAPTER  IV 
RESULTS 

The  specific  computer  program  which  resulted  from  this  study  is 
found  in  Appendix  B.    As  each  portion  of  the  program  was  developed,  it 
was  tested  in  such  a  way  that  every  part  of  the  program  was  operated. 
Intermediate  print-outs  of  various  steps  were  obtained  and  verified  for 
all  possible  types  of  inputs  for  that  particular  portion.,     Efforts  were 
directed  toward  simplifying  and  shortening  the  routines  involved. 

As  the  final  program  was  amalgamated,  the  testing  and  review 
continued.    When  it  was  found  to  be  functioning  properly,    the  various 
test  print-out  statements  and  corresponding  comment  cards  were  removed. 
Only  those  comment  cards  which  will  assist  in  reviewing  the  program  remain. 

Sample  illustrations  of  a  distress  case  involving  an  aircraft  (ditching) 
with  associated  survival  craft  ranging  from  life  jacket,,  life  raft  to  life 
boat,  with  a  ship  as  the  search  unit  are  found  in  the  following  section. 
These  problems  show  the  input  data,  the  manual  solution,  and  the 
computer  solution.     Manual  solutions  were  done  on  small  area  plotting 
sheets.     It  was  noted  that  when  the  manual  solutions  were  done  by  people 
thoroughly  familiar  with  the  problem,   human  error  lead  to  one  17  mile 
error  in  datum.    The  computer  solution  was  correct. 

The  running  time  for  the  samples  varied  with  the  number  of  the  search. 
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The  longest  time,   68  seconds,  occurred  on  the  run  for  the  initial  search, 
and  54  seconds  were  required  for  the  subsequent  searches.     The  first 
search  has  a  longer  print-out  associated  with  the  distress  message. 
It  also  involved  computations  of  initial  position  error  of  the  distress 
incident  and  SAR  unit  navigation  error  as  well  as  the  immersion  hypothermia 
and  chill  factor  relationships  which  were  not  computed  in  subsequent  runs. 
All  r/fmning  times  would  be  reduced  if  this  program  were  in  an  operational 
status  as  the  various  tables,   program  instructions  and  other  associated 
data  would  already  be  in  storage  within  the  computer.     This  information 
had  to  be  entered  at  the  beginning  of  each  of  these  sample  runs. 
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I.      TRIAL  DATA  COMPUTATIONS 


RUN  ONE 


INPUT  DATA 

Classification  of  distress  unit  position  report 

Classification  of  SAR  unit  navigation  method 

Initial  reported  distress  position 

Position  of  distress  unit's  last  fix 

SAR  unit  position  at  start  of  case 

Time  of  distress  incident 

Time  specified  for  arrival  of  SAR  unit 

Surface  wind  at  distress  position 

Surface  wind  for  past  24  hours  at  distress 

position 
Average  surface  current  at  distress  position 
Air  temperature  at  distress  position 
Sea  water  temperature  at  distress  position 
Search  number 
Type  lifesaving  craft  employed 


Aircraft  -  DR 
Ship  -  FIX 
42.  50N,    52.42W 
41.50N,    61.50W 
50.00N,   35.00W 
301730Z  (JAN) 
021130Z  (FEB) 
340°T  @  17  lets. 

265°T  @  10  kts0 

195°T@  15.6  miles/day 

12.8°C 

7.2°C 

1 

Life  raft 
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RESULTS  BY  MANUAL  COMPUTATIONS 


Distress  unit  navigation  error  (X) 
SAR  unit  navigation  error  (Y) 
Forecast  datum  of  survivors 
Forecast  search  area  limits 


51.3  miles 

5.  0  miles 

40.97N,    51.80W 

42. 14N,    39.81N 
50.25W,    53.35W 


RESULTS  BY  COMPUTER  (PROGRAM)  COMPUTATIONS 


Distress  unit  navigation  error  (X) 
SAR  unit  navigation  error  (Y) 
Forecast  datum  of  survivors 
Forecast  search  area  limits 


51.1  miles 

5.  0  miles 

40o90N,    51.84W 

42.01N,    39.80N 
50.37W,    53.30W 


INPUT  DATA 


RUN  TWO 


Last  computed  distress  position 

Type  lifesaving  craft  employed 

Time  associated  with  last  distress  position 

Time  specified  for  next  search  to  begin 

Surface  wind  at  computed  distress  position 


40.90N,    51.84W 
Man  in  water 
021130Z  (FEB) 
041330Z  (FEB) 
180°T  @  20  kts. 
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Surface  wind  for  past  24  hours  at  computed 

distress  position 
Average  surface  current  at  computed  distress 

position 
Search  number 

RESULTS  BY  MANUAL  COMPUTATIONS 


340°T  @  17  kts. 


195°T  @  15.  6  miles/day 


Forecast  datum  of  survivors 


Forecast  search  area  limits 


39.88N,    52.  16W 

41.30N,    38.46N 
50.28W,    54.04W 


RESULTS  BY  COMPUTER  (PROGRAM)  COMPUTATIONS 


Forecast  datum  of  survivors 


Forecast  search  area  limits 


39.88N,    52.14W 

41.30N,    38.47N 
50.30W,    53.98W 


RUN  THREE 


INPUT  DATA 


Last  computed  distress  position 

Type  lifesaving  craft  employed 

Time  associated  with  last  distress  position 

Time  specified  for  next  search  to  begin 


39.88N,   52.  14W 
Life  boat 
041330Z  (FEB) 
062130Z  (FEB) 
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Surface  wind  at  computed  distress  position 
Surface  wind  for  past  24  hours  at  computed 

distress  position 
Average  surface  current  at  computed  distress 

position 
Search  number 


095OT  @  04  kts. 

180°T  @  20  kts. 

180°T  @  15.  6  miles/day 
3 


RESULTS  BY  MANUAL  COMPUTATIONS 


Forecast  datum  of  survivors 


Forecast  search  area  limits 


39.89N,   52.09W 

41.62N,    38.16N 
49.83W,    54.35W 


RESULTS  BY  COMPUTER  (PROGRAM)  COMPUTATIONS 


Forecast  datum  of  survivors 


Forecast  search  area  limits 


39.87N,    52.07W 

41.60N,    38.14N 
49.82W,    54.32W 


SUMMARY  OF  RUNNING  TIMES 


Run  number  one 


Run  number  two 


Run  number  three 


1  min.   8  sec. 


0  min.   54  sec, 


0  min.   54  sec. 
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CHAPTER  V 
SUMMARY  AND  CONCLUSIONS 

L       SUMMARY 

Growing  maritime  operations  have  found  the  need  for  more  effective 
search  and  rescue  facilities  expanding.     One  area  of  search  and  rescue 
operations,   that  of  the  determination  of  search  areas  and  datum  points, 
has  been  studied  to  ascertain  the  feasibility  of  the  use  of  computers  in 
solving  this  problem.     Using  the  National  Search  and  Rescue  Manual  as 
a  guide  to  the  present  doctrine  employed  by  this  country,   a  computer 
program  was  developed,  tested  and  found  to  be  an  operationally  feasible 
program. 

During  the  course  of  the  development  of  this  program  certain  areas 
of  the  current  doctrine  Were  reviewed  in  the  light  of  recent  literature  and 
theories  which  may  offer  improved  techniques.     The  conclusions  reached 
in  this  review  are  discussed  in  the  following  section. 

II.      CONCLUSIONS 

The  basic  doctrine  outlined  in  the  National  Search  and  Rescue  Manual 
is  amenable  to  computer  programming.    Although  the  program  which  was 
developed  does  not  cover  every  possible  aspect  of  this  doctrine,   it  does 
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handle  one  of  the  basic  types  of  search  problems  frequently  encountered. 
This  program  can  greatly  reduce  the  time  required  to  solve  the  problem 
associated  with  the  location  of  datum  and  its  surrounding  search  area. 
The  accuracy  of  the  computer  solution  is  dependent  entirely  upon  the 
accuracy  of  the  input  information.     It  greatly  reduces  the  possibilities  of     i 
human  error  in  carrying  out  the  rather  extensive  calculations  involved  in 
the  problem0    As  new  methods  are  developed  for  determining  various 
portions  of  the  input  data  with  greater  accuracy,  the  value  of  the  computer 
solution  will  greatly  increase. 

In  working  with  the  National  Search  and  Rescue  Manual  certain 
conclusions  were  reached.     The  first  of  these  was  that  the  doctrine 
.outlined,,  although  not  the  ultimate  in  sophistication,   certainly  represented 
an  excellent  model  for  use  in  existing  RCC  facilities  which  were  not 
equipped  with  automatic  data  processing  equipment  or  with  large  staffs. 
Admittedly  there  were  areas  of  weakness  in  the  basic  doctrine.     However, 
to  refine  the  doctrine  to  fully  consider  the  different  shortcomings  would 
have  unnecessarily  complicated  the  problems  at  the  sacrifice  of  valuable 
time.    And  there  is  no  definite  measure  of  how  this  additional  refinment 
would  alter  the  risks  associated  with  the  resultant  decisions. 

However,  the  review  of  the  literature  in  this  study  did  suggest 
several  possible  areas  for  improvement  if  a  computer  installation  is 
available,  and  one  suggestion  is  applicable  with  existing  manual  operations, 
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Considering  the  computer  applications  first,   it  is  apparent  that  the 
problem  of  free  fall  and  parachute  drift  is  subject  to  operational  analysis 
by  computers.     The  present  running  time  of  slightly  less  than  ten  minutes 
associated  with  Cohan's  program  can  probably  be  reduced  as  his  program 
is  refinedo 

Also,  a  computer  could  use  the  program  developed  by  Burns  at  the 
Navy  Oceanographic  Office  to  determine  a  drift  elipse  area  based  upon 
the  statistical  current  roses  found  in  the  Atlas  of  Surface  Currents. 

The  development  by  the  Fleet  Numerical  Weather  Facility  of  the 
system  for  forecasting  surface  currents  appears  to  have  the  greatest 
potential  for  improving  present  doctrine  in  that  it  may  be  used  either  with 
or  without  a  computer.     The  program  developed  by  this  study  can  be 
modified  to  receive  information  from  the  Fleet  Numerical  Weather  Facility 
on  the  predicted  surface  currents  at  the  four  points  of  the  grid  containing 
datum.     These  vectors,  which  reflect  a  combination  of  the  permanent 
circulation  pattern  and  local  wind  currents,   can  then  be  resolved  into  a 
vector  for  the  datum  point.     Since  this  information  is  based  on  recently 
reported  weather  conditions  and  sea  surface  temperatures,  it  should  tend 
to  be  more  accurate  than  the  information  taken  from  the  Atlas  of  Surface 
Currents  which  represents  a  statistical  summary  of  many  years  based,  in 
s.ome  cases,  on  a  single  observation. 

This  analysis  of  Fleet  Numerical  Weather  Facility  information  cquld 
also  be  graphically  determined  or  manually  calculated  without  the  assistance 
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of  a  computer0 

Laevastu's  work  indicates  that  local  wind  current  divergence  can  be 
treated  in  a  more  refined  manner  by  considering  small  increments  of 
latitude  and  wind  velocities.     This  refinement  could  be  handled  very 
easily  by  modifying  that  portion  of  the  specific  program  relating  to  this 
computation, 

III.       IMPLICATIONS 

The  program  developed  in  this  study  can  be  rewritten  in  the  computer 
language  used  in  the  AMVER  Center,   Office  of  the  Commander  Eastern 
Area,     Thus,   it  can  serve  to  assist  the  handling  of  these  distress  cases 
which  occur  in  the  N0rth  Atlantic.     It  is  suggested  that  serious  consideration 
be  given  to  the  adaptation  of  the  Fleet  Numerical  Weather  Facility  data 
services  and  the  associated  modification  of  this  basic  program.     Should 
this  latter  suggestion  not  be  practical,   consideration  should  be  directed 
toward  utilizing  Burns'  program  in  conjunction  with  the  Atlas  of  Surface 
Currents,    As  the  Pacific  AMVER  system  becomes  fully  operational,  this 
this  program  can  be  further  modified  to  handle  SAR  cases  falling  under  the 
jurisdiction  of  Commander  Western  Area. 

In  view  of  the  number  of  entries  found  in  the  Atlas  of  Surface  Currents 
(approximately  20,  000  for  the  Atlantic  and  40,  000  for  the  Pacific)  it  is  not 
recommended  that  this  information  be  placed  in  storage  and  the  basic 
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program  modified  to  handle  that  portion  of  the  Average  Sea  Current  problem 
automatically.     Individual  cases  should  continue  to  use  this  information  as 
manually  extracted  from  the  appropriate  charts. 

IV.      RECOMMENDATION  OF  AREAS  FOR  FURTHER  STUDY 

This  study  has  served  as  a  first  step  in  the  development  of  a  full-scale 
computerized  SAR  program.    As  such,   it  has  opened  up  many  new  areas  where 
further  study  and  research  should  prove  fruitful.     Some  of  these  areas  are 
discussed  briefly  below. 

Determination  of  a  datum  and  associated  search  area  when  the  time  of  the 
incident  and  estimated  position  are  not  known. 

When  a  ship  or  aircraft  is  missing  or  overdue,   the  RCC  Controller 
is  faced  with  the  problem  of  estimating  the  probable  boundaries  of  an 
area  in  which  the  distress  may  have  occurred.     He  must  establish  the 
limits  of  the  time  period  in  which  the  distress  could  have  occurred.    With 
this  information  he;  then,  applies  minimum  and  maximum  possible  values  of 
drift  for  the  time  period  in  question  starting  with  the  earliest  possible 
distress  time  and  ending  with  the  latest  possible  distress  time. 

Development  of  a  fully  operational  free  fall  and  parachute  drift  program. 

Refinement  and  further  development  of  the  work  done  by  Armacost, 

80 


Saunders  and  Cohan  can  materially  contribute  to  the  effectiveness  of  a 
computer  program  for  SAR  operations. 

Development  of  a  program  to  handle  small  craft  (inshore)  distress  cases. 

With  the  increasing  interest  in  recreational  boating  displayed  by 
the  general  public,  it  is  mandatory  that  improved  techniques  be  developed 
to  handle  the  associated  increase  in  SAR  cases  occurring  with  these  craft. 
It  has  been  estimated  that  approximately  85%  of  all  SAR  cases  handled  by 
the  U.  S.  Coast  Guard  occur  within  twenty  to  twenty-five  miles  of  the 
coast.    Someplace  in  the  neighborhood  of  95%  of  all  cases  are  found  to 
take  place  within  100  miles  of  the  coastline.  48 

One  of  the  major  difficulties  in  treating  these  cases  is  associated 
with  the  effects  of  tidal  currents ,  rotary  currents,  and  run-off  currents , 
all  of  which  are  effected  by  the  configuration  of  the  coastline. 

Development  of  a  program  relating  search  area,   search  units,  and  search 
patterns . 

Once  a  search  area  has  been  specified,  the  RCC  Controller  is  faced 
with  the  problem  of  recommending  search  patterns  for  those  available 


48 

Bo  R.  House.  Analysis  of  U.  S.  Coast  Guard  Aviation  Operations 

(Report  Number  SER-50371,  Stratford:  Sikorsky  Aircraft  Division  of  United 

Aircraft  Corporation,    1964),   p.  36. 
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search  units  which  are  participating  in  the  case.    Since  the  type  of  search 
pattern  recommended  is  a  function  of  the  number  and  type  of  search  units 
participating,  the  weather  and  other  operating  conditions  at  the  scene, 
and  the  various  probabilities  of  detection  desired,  the  development  of  a 
program  to  assist  in  making  this  recommendation  would  be  most  beneficial. 
In  essence  this  would  be  a  program  treating  the  doctrine  outlined  in 
Chapter  Seven  of  the  National  Search  and  Rescue  Manual. 

Studies  evaluating  the  frequency  of  distress  cases  and  the  location  of  SAR 
facilitieso 

As  the  new  system  of  reports  for  distress  cases  becomes  fully 
operational,  actual  historical  data  will  become  available.    Analysis  of 
this  data  and  existing  SAR  facilities  should  shed  new  light  on  possible 
relocation  and  expansion  of  SAR  facilities.     Studies  similar  to  that  performed 
by  DuPeza,   Eustis,  Guthrie  and  Zook^^  should  be  of  great  value. 

Continued  drift  analysis. 

At  the  present  time  the  U.  S.  Coast  Guard  has  a  drift  datum  marker 
bfctoy  which  is  being  evaluated  by  operational  field  testing.     The  data 
obtained  from  these  tests  coupled  with  the  current  studies  being  conducted 

49 

Jules  B.   DuPeza  and  others,    "A  Method  for  Evaluating  SAR  Requirements 

for  the  Twelfth  Coast  Guard  District"  (Unpublished  Master's  research  paper, 

The  United  States  Naval  Postgraduate  School,   Monterey,    191 
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by  the  U.  S.   Naval  Oceanographic  Office,   U.   S.   Naval  Fleet  Numerical 
Weather  Facility,   and  other  oceanographers  should  provide  a  great  deal 
of  information  which  can  lead  to  a  better  understanding  of  the  problems 
associated  with  ocean  current  analysis. 

In  all  cases,  the  continued  application  of  man's  growing  knowledge 
and  utilization  of  the  tools  which  he  has  developed  can  be  directed 
toward  the  benefit  of  mankind.     The  success  or  failure  of  such  efforts 
is  limited  only  by  the  imagination,   enthusiasm  and  patience  of  those 
who  truly  seek  the  answers. 
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APPENDIX  A. 


FLOWCHART  SYMBOLS  AND  DEFINITIONS 


FLOWCHART  SYMBOLS  AND  DEFINITIONS 


-     Terminal  indicator 


-     Input/output  indicator 


Decision  indicator 


Processing,  annotating  indicator 


0 


-     Off  page  connector 


On  page  connector 


O   O    V  A  ~     plow  indicator 
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APPENDIX    B. 


PROGRAM  SAR 


Program  SAR  is  the  specific  program  developed  in  this  study.     It 
determines  the  datum  point  and  associated  search  areas  for  distress 
cases  involving  ships  and  aircraft  which  ditch.     This  Appendix  includes 
the  definitions  of  the  terms  and  inputs  found  in  the  program,   the 
associated  computer  program  itself,  and  sample  print-out  as  discussed 
at  the  end  of  Chapter  III. 
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DEFINITIONS  OF  ALL  INPUT  VARIABLES 


VARIABLE 
NAME 


ATEM 


CRAFT 


DCLS 


DDATE 

DFTERR 

DHR 

DLAT 

DLONG 

DMIN 

DMO 


DEFINITION 

Air  temperature  at  distressed  position  in  degrees  centigrade 
(°C). 

Type  lifesaving  craft  that  RCG  Controller  thinks  survivors 
will  employ: 

Life  jacket  -  code    1. 

Boat  -  code    2^, 

Raft  -  code   _3._ 

Classification  of  Distressed  unit's  position  report: 

Surface  DR  -  code  01 

Aircraft  DR  -  code  02 

Surface  FK  =  code  .05. 

Aircraft  FDC  -  code  J_0 
Surface  and  Aircraft  DF 

Class  A  -  code  05. 

Class  B  -  code  20 

Class  C  -  code  50 
Surface  and  Aircraft  DF-FCC 

Class  A  -  code  20. 

Class  B  -  code  40 

Class  C  -  code  60 

Class  D  -  code  80 

Date  of  distress  in  ZULU  (GMT)  time. 

Drift  error  from  Min-Max  Table  in  CG-308. 

Hour  of  distress  in  ZULU  (GMT)  time. 

Latitude  of  distress  in  degrees  and  hundredths. 

Longitude  of  distress  in  degrees  and  hundredths. 

Minute  of  distress  in  ZULU  (GMT)  time. 

Number  of  days  in  the  month  associated  with  DDATE. 
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DROG  Drift  of  a  raft  or  boat  without  a  drogue  as  a  factor  of  wind 

velocity.     Values  obtained  from  a  graph  in  CG-308. 

EXH  Time  until  an  immersed  man  will  become  exhausted  or 

unconscious  as  a  factor  of  sea  water  temperature.     Obtained 
from  table  in  Atlas  of  Surface  Currents. 

FLAT  Latitude  of  last  fix  of  distressed  unit  in  degrees  and 

hundredths. 

M  Number  of  the  search  that  is  being  computed. 

PWDIR  Wind  direction  at  the  distress  position  for  the  24  hour  period 

prior  to  the  distress. 

PWIND  Wind  velocity  at  the  distress  position  for  the  24  hour  period 

prior  to  the  distress „ 

SCLS  Classification  of  the  SAR  unit's  position  determination: 

Surface  DR  -  code  03. 

Aircraft  DR  -  code  04. 

Surface  FIX  -  code  05 

Aircraft  FLX  -  code  10 
Surface  and  Aircraft  DF 

Class  A  -  code  05. 

Class  B  -  code  20. 

Class  C  -  code  50 
Surface  and  Aircraft  DF-FCC 

Class  A  -  code  20. 

Class  B  -  code  40 

Class  C  -  code  60. 

Class  D  -  code  80 

SCURD  Direction  of  surface  current  at  distress  position  taken  from 

the  Atlas  of  Surface  Currents. 

SCURV  Velocity  of  surface  current  at  distress  position  taken  from 

the  Atlas  of  Surface  Currents . 

SDATE  Date  specified  by  RCC  Controller,   usually  the  ETA  of  the  SAR 

unit  on  scene  in  ZULU  (GMT)  time. 

SEARCH         Safety  factor  associated  with  'M'. 
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SHR  Hour  specified  by  RCC  Controller,  usually  the  ETA  of  the  SAR 

unit  on  scene  in  ZULU  (GMT)  time. 

SLAT  Latitude  of  the  SAR  unit's  last  fix  in  degrees  and  hundredths. 

SLOW  Longitude  of  the  SAR  unit's  last  fix  in  degrees  and  hundredths. 

SMIN  Minute  specified  by  the  RCC  Controller,  usually  the  ETA  of 

the  SAR  unit  on  scene. 

STEMP  Temperature  of  the  seawater  at  the  distress  position^  °C. 

STM  Survival  time  for  a  man  immersed  in  water  at  a  certain 

temperature  as  a  factor  of  the  temperature. 

WDIR  Wind  direction  at  distress  position  in  degrees  True. 

WIND  Velocity  of  wind  at  distress  position. 

WINTEM  Wind  chill  factor  for  the  effect  of  wind  on  exposed  flesh. 

X  Position  error  for  the  distress  report. 

Y  Position  error  for  the  SAR  unit. 
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INPUT  DATA  FOR  FIRST  SEARCH 


VARIABLE 

DATA  CARD 

NAME 

COLUMNS 

DCLS 

1,2 

SCLS 

3,4 

SOURCE  AND  FORMAT 
OF  INPUT  DATA 


DLAT 


DLONG 


FLAT 


FLON 


5-10 


11  =  16 


17-22 


23-28 


Classification  of  distressed  unit's  position 
report.     Two  digit  code  number. 

Classification  of  search  unit's  position  report. 
Two  digit  code  number. 

Latitude  of  distress  incident.     Six  digit  number 
of  degrees  and  hundredths  (i.e0   45°39'  =  045.65). 

Longitude  of  distress  incident.     Six  digit  number 
of  degrees  and  hundredths  (Le.   75°21  -  075.  35). 

Latitude  of  last  fix  of  distressed  unit  six  digit 
number  of  degrees  and  hundredths,     (i.e. 
42°54'=042.90). 

Longitude  of  last  fix  of  distressed  unit.  Six 
digit  number  of  degrees  and  hundredths  (i.  e. 
63°13'=063.22). 


SLAT 


SLON 


CRAFT 


DDATE 


DHR 


29-34  Latitude  of  SAR  unit  at  start  of  incident.    Six 

digit  number  in  degrees  and  hundredths  (i.e. 
40°38"=040.63). 

35-40  Longitude  of  SAR  unit  at  start  of  incident.    Six 

digit  number  in  degrees  and  hundredths  (i.  e. 
68°07'=068.12). 

41,  42  Type  of  life  saving  device  that  the  RCC  coordinator 

believes  survivors  will  be  using.     One  digit  code 
number  followed  by  a  decimal. 

43,  44  Date  of  distress  incident.    Two  digit  number  from 

the  distress  DTG. 

45,  46  Hour  of  the  distress  incident.     Two  digit  number 

from  the  distress  DTG. 
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VARIABLE 

DATA  CARD 

NAME 

COLUMNS 

DMIN 

47,48 

SDATE 

49,50 

SOURCE  AND  FORMAT 


OF  INPUT  DATA 


SHR 


SMIN 


WIND 


PWIND 


WDIR 


PWDIR 


SCURV 


SCURD 


51,52 

53,54 

55,56 
57,58 
59-61 

62-64 

65=68 


69-71 


Minute  of  the  distress  incident, 
number  from  the  distress  DTG. 


Two  digit 


Date  specified  by  the  RCC  Controller,  usually 
the  date  (GMT)  that  the  SAR  unit  will  arrive  on 
scene.     Two  digit  number. 

Hour  specified  by  the  RCC  Controller,  usually 
the  hour  (GMT)  that  the  SAR  unit  will  arrive  on 
scene.     Two  digit  number. 

Minute  specified  by  the  RCC  Controller,  usually 
the  minute  that  the  SAR  unit  will  arrive  on  scene. 
Two  digit  number. 

Wind  velocity  at  distress  position  from  the  latest 
weather  information.     Two  digit  number. 

Wind  velocity  for  the  previous  24  hours  at  the 
distress  position.     Two  digit  number. 

Direction  from  which  the  wind  is  blowing  at  the 
distress  position  at  time  of  distress.     Three 
digit  number. 

Direction  from  which  the  wind  had  been  blowing 
at  the  distress  position  for  the  24  hour  period 
prior  to  the  distress  time.     Three  digit  number. 

Surface  current  velocity  at  the  distress  position 
taken  from  the  Atlas  of  Surface  Currents.    Four 
digit  number  (nautical  miles  per  day)  in  miles 
and  tenths  (i.e.   09.9). 

Direction  of  the  surface  current  flow  at  the 
distress  position  taken  from  the  Atlas  of  Surface 
Currents.     Three  digit  number. 
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VARIABLE        DATA  CARD 


NAME 


COLUMNS 


SOURCE  AND  FORMAT 
OF  INPUT  DATA 


ATEM 


72-75  Air  temperature  at  the  distressed  position  at  the 

distress  time  taken  from  weather  reports  in  degrees 
centigrade.     Four  digit  number  in  degrees  and 
tenths  (i.  e.   17.  6). 


STEMP 


76-79  Sea  water  temperature  at  distress  position  in 

degrees  centigrade  taken  from  latest  weather 
data.     Four  digit  number  in  degrees  and  tenths 
(i.e.    15.3). 


M 


80 


Search  number  for  which  datum  is  to  be  computed. 
One  digit  number. 


DMO 


1,2 


Number  of  days  in  the  month  of  DDATE.     Two 
digit  number. 


X 

Y 


3-7 
8-12 


Left  blank  for  the  first  search. 


Left  blank  for  first  search. 
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INPUT  DATA  FOR  SECOND  THRU  FIFTH  SEARCHES 


VARIABLE        DATA  CARD 


NAME 
DCLS 
SCLS 
DLAT 

DLONG 

FLAT 

FLON 

SLAT 

SLON 

CRAFT 

DDATE 

DHR 

DMIN 

SDATE 

SHR 

SMIN 


COLUMNS 
1,2 
3,4 
5-10 

11  =  16 

17-22 
23-28 
29-34 
35-40 
41,42 
43,44 
45-46 
47,48 
49,50 

51,52 

53,54 


SOURCE  AND  FORMAT 
OF  INPUT  DATA 


Blank. 
Blank. 

DTMLAT  from  preceding  search  results.    Six 
digit  number  in  same  format  as  DTMLAT. 

DTMLONG  from  preceding  search  results.    Six 
digit  number  in  same  format  as  DTMLONG. 

Blank. 

Blank. 

Blank. 

Blank. 

Same  as  CRAFT  in  preceding  search. 

SDATE  from  preceding  search. 

SHR  from  preceding  search. 

SMIN  from  preceding  search. 

New  specified  date  by  RCC  Controller,  usually 
the  estimated  date  (GMT)  that  the  SAR  unit  will 
complete  the  preceding  search. 

New  specified  hour  by  the  RCC  Controller,  usually 
the  estimated  hour  (GMT)  that  the  SAR  unit  will 
complete  the  preceding  search. 

New  specified  minute  by  the  RCC  Controller,  usually 
the  estimated  minute  (GMT) that  the  SAR  unit  will 
complete  the  preceding  search. 
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VARIABLE        DATA  CARD 


NAME 


COLUMNS 


SOURCE  AND  FORMAT 
OF  INPUT  DATA 


WIND 


PWIND 


55,  56  Wind  velocity  at  the  preceding  computed  datum 

position.,     Taken  from  weather  reports  for  the 
preceding  specified  time  (SDATE,   SHR,   SMIN). 
Two  digits  in  same  format  as  search  one. 

57,  58  Wind  velocity  for  the  previous  24  hours  period 

at  datum  computed  in  the  preceding  run.     Two 
digit  number  obtained  from  weather  reports. 


WDIR 


PWDIR 


SCURV 


59-61  Wind  direction  at  datum  computed  in  the  preceding 

run.     Three  digit  number  obtained  from  the  latest 
weather  report. 

62-64  Direction  of  previous  24  hour  wind  at  datum 

computed  in  preceding  run.     Three  digit  number 
obtained  from  weather  reports. 

65-68  Surface  current  velocity  at  datum  computed  in 

preceding  run.     Four  digit  number  of  miles  and 
tenths  (i.  e.   13.  4)  obtained  from  the  Atlas  of 
Surface  Currents. 


SCURD 


69-71  Direction  of  surface  current  at  datum  computed  in 

preceding  run.     Three  digit  number  obtained  from 
Atlas  of  Surface  Currents. 


ATEM 
STEMP 
M 
DMO 


72-75  Blank. 

76-79  Blank. 

80  Number  of  the  search,   2  thru  5.     One  digit. 

1,  2  Number  of  days  in  the  month  of  DDATE.    Two 

digit  number. 

3-7  Value  of  DX'  from  the  initial  search.    Five  digit 

number  with  one  decimal  place  o 

8-12  Value  of'Y'  from  the  initial  search.     Five  digit 

number  with  one  decimal  place. 
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FIRST  SEARCH  PRINT  OUT 
XXXXXXX        Z 
FM    COMEASTAREA 
TO    XXXXXXXXXXXXXXXXXXXXX 
BT 

UNCLASS 
A.     DISTRESSED  XXXXXXXXXXXXXXX 

1.  FORECAST  DATUM  AT  02113QZ.     40.90N,     51.84W. 

2.  FORECAST  SEARCH  AREA  BOUNDED  BY    42.  01N„     39.8QN. 
50.37W,     53.30W. 

3.  SURVIVOR  CONDITION  FROM  30173QZ  OF  MAN  IN  WATER. 
TIME  UNTIL  EXHAUSTION      1.  00  HOURS.     TIME  UNTIL 
DEATH     3.00  HOURS. 

4o     EFFECT  OF  SURFACE  WIND  SAME  AS    -27.0  DEGREES 

FAHRENHEIT  WITH  ZERO  WIND. 
5o 


BT 

RELEASED  BY  XXXXXXXXXXXXXXXXXXXX 
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SUBSEQUENT  SEARCH  PRINT  OUT 

FM    COMEASTAREA 

TO    )QQQQQQCXXXXXXXXXXXXXX 

BT 

UNCLASS 

A,     DISTRESSED  XXXXXXXXX3QCXXXXXX 

1.  FORECAST  DATUM  AT  082130Z.      39.87N,     52.  07W. 

2.  FORECAST  SEARCH  AREA  BOUNDED  BY    41.  6QN,     33.  14^ 
49.82W,     54. 32W. 


BT 


RELEASED  BY  XXXXXXXXXXXXXXXXXXXX 
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APPENDIX  C 

"Parachute  Drift" 

This  Appendix  is  an  extract  from  the  study,  "Computer 
Program  of  USCG  Search  and  Rescue  Procedures,"  which  was  an 
unpublished  undergraduate  study  performed  by  Cadets  First 
Class  Robert  L.  Armacost  and  Norman  T.  Saunders  at  the  United 
States  Coast  Guard  Academy,  New  London,  Connecticut,  in  May 
1964. 
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RLA  PARACHUTE  DRIFT 

10  Feb.,  1964 

After  examining  the  Parachute  Drift  Table,  Figure  6-3 
in  GG-303,  it  seemed  that  most  of  the  tabular  values  were 
linearly  related.   I  plotted  these  values  as  a  family  of 
curve 3  of  wind  speed.  I  used  the  height  of  the  parachute 
opening  as  the  ordinate,  and  the  actual  drift  as  the  abscissa. 
The  curves  didn't  prove  to  be  exactly  linear  but  they  were 
very  nearly  so.   For  some  reason,  I  don't  feel  that  this 
linearity  should  be  so  apparent.   I  decided  to  try  to  arrive 
at  some  purely  analytical ; equation  to  describe  the  motion  of 
the  man  as  he  comes  down. 

If  we  assume  that  he  is  jumping  from  a  great  height,  we 
know  that  he  has  to  free  fall  to  at  least  IS, 000  feet.   At 
higher  altitudes,  he  may  suffer  from  anoxia,  or  stands  the 
chance  of  freezing  since  the  temperature  is  -5.3°F.   In  order 
to  consider  his  motion  in  free  fall,  we  must  examine  the  com- 
ponents of  his  motion.   They  are  due  tor 

1.  The  force  of  gravity. 

2.  The  effect  of  the  wind  on  him  in  a  horizontal  - 
direction. 

3.  His  motion  due  to  the  velocity  of  the  airplane. 
All  of  these  components  are  opposed  by  the  air  resist- 
ance. The  need  for  these  equations  is  so  we  may  have  a  very 
accurate  entry  for  the  computer. 

The  derivation  follows. 
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Vertical  Component  of  Free  Fall 

Let  the  y-axis  be  the  vertical  axis  and  be  positive  down. 

Air  Resistance,  R  -  c  v^. 

Definitions: 

ra  =  mass  of  the  man 

g  9   acceleration  due  to  gravity  (assumed  constant) 

=32.2  ft/sec2 

W  -  weight  of  the  man  s  mg 

am  s  dvm  s  acceleration  of  the  man 
"HE"" 

R  5  air  resistance 

c  s  factor  of  proportionality 

We  may  now  start  to  write  our  equations. 

radvjx  »  mg  -  cr2^2 

Gt 

Let  k2  =•  fi2,,  then  l.dvja.  *  g^  -  vm^ 
m  It  ^"^    Ic 

Let  u2  Z  g/k2 

2 

dvm   a  k  dt 

~H — Z 

vQ   =  0  at  time  t  »  0 


dvm   n  f  k^dt 


J  ur-  vm 
yo 

_l.tanh  (vm/u)  =  k  t  -^.t 
u  -  u2 


vm  *  u»tanhC£*t)  (l) 
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As  t">°«,  vm  s  u  as  a  limiting  velocity. 
We  may  write  further  equations  of  motion: 

dv  -  u>tanh(§*t) 

dt 

At  t  =  0,     y  *  0 


i 


dv  -  [  u»tanh(  g  »t) 
dt  J  u 


y  s  u  In  coshl-S— *t} 

y  -  u^ln  cosh(-J-*t) 
g 

u 

Then:  y  s  u  In  cosh(Xt) 
A 

A  y  s  In  cosh  (At) 
u 


W 


=  coshCxt) 
t  s  l„cosh   e 


(2) 


(3) 


Refer  to  page  51,  Betz,  Burnham,  Ewing.   They  say  u  -   175ft. 

sec. 

for  a  man  in  free  fall.   From  this  we  may  evaluate  an  approxi- 
mate value  of  c. 


u  = 


jg. 


k  z  _£_ 


m 
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u  3  Vgm     -   ywm 

*  c 


u  175 

c  =  0.0766  (A) 

Horizontal  Component  of  Motion  Due  to  Wind 

Assume:  1.  Horizontal,  steady  wind 

2.   The  average  man  is  6  feet  tall  and  weighs  ISO 
pounds.  He  presents  a  projected  area,  k^   -  1125  in? 

Neglecting  air  resistance,  we  can  write  a  few  funda- 
mental equations  and  find  an  equation  for  this  motion, 

F  ~  force  (pounds) 

p  s  pressure  (pounds  per  square  inch) 

AjQ  s  projected  area  of  the  man  (ft  ) 

V^  3  horizontal  velocity  of  the  man 

Vw  r  velocity  of  the  wind 

F  =  P^m  *  mah  «  nidVjj. 

dt 

P  =■  y^r    where  velocity  head  3  V~. 
zg  zg 

and  density  (pounds/ft  J  is  equal  to  Y 

%i4^k  =  XVw  Am 
dt     zg 

But  1%  =  Wm 

•    dVm   ~  ^iZ.Am  -    ^Yw2Am 
dt      zgWrn./        zWm 
g 
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Let  K2  -  y^Am 

zWm 

Then:      K  s  Vw^/xK£Z     =  Vi/0«0gU^3F 

K  »   (4.16  x  1(T2)VW  (B) 

K2  s  17.33  x  10-HW2  (C) 

dVh   =  K2 
dt 

Vh  -  K2t  .  17.33  x  lCT4Vw2t 

For  a  wind  velocity  of  60  knots   (101  ft/sec.): 

Vh  =  (17.33  x  10-^1 (I0l)2t 

For:        V  <   15  ft/sec,  R  =.  cv 

15  ft/sec  <  V  <  Sonic,   R  =  c2^2 
To  get  some  idea  of  our  proportionality,  let  us  assume  a 
free  fall  of  2000  feet;  we  must  compute  the  time  and  thus,  V^. 

y  s  £000 

c  s  0.0766     c2  s  53. S  x  10"4 

m-  =  W^  3  130'   s  5.59 
g    32.2 

k2  s  c2/m  =10.5  xlO"^ 

u2  s  g/k2  =  3.06  x  IcA 

X  =  g/u  *  32.2/175  =  0.184 

Substituting  these  values  back  into  Equation  (3)  and  solving: 

.  t  =  15.2  seconds 
and    Vh  -  17. 7t  =269  ft/sec.  =  159.5  knots 
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Therefore,  we  may  assume  that  the  resistance  to  his 

horizontal  motion  is  proportional  to  th©  square  of  the 

velocity.  Since  he  is  in  the  same  medium,  the  constant  of 

proportionality  should  be  the  same  since  he  should  have  the 

same  limiting  velocity  if  sufficient  force  were  applied. 

dVh  -  K2-c2Vh2 
dt 

" . :  ~£k '.     stfb 

K2-  c2Vh2 

„  dVh  =  c2dt 

K2  -  V  2 
c2 

Let  B2  -  K2/c2 

Then:     d^h   m   c2dt 
B2  -  Vh2 

At  t  z  0,  %  ~   0 

r  Vh  o  -   ftc2dt 
J0  B^Vh2       Jo 

1  tann^C  ¥h/B)   -  c2t  «  K2t 
B  IF 

Vh  r  B.tanh(K2t/B)  (4) 

As  t  increases,  Vn  approaches  B  as  a  limiting  velocity. 

B  = JL.  =  4rX6  x  lQ-% 
c         7.66  x  10"^ 

B  =  0.543VW 

-x 


f-&=    |    Btanh(K2t/B) 
Jodt       Jo 
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At  t  =  0,   x 


x   =  B2ln  cosh(K2t/B) 


~F 


But:    B2  -  1 

F    "? 


x  =     LJin  coshlK^/B)    =     1  In   cosh /l7. 33  x  KT^vSt^ 
"c^  c2  \0.543VW  / 

X  -     1  In  cosh(3.19  x  10"3Vwt)  (5) 

c 

CG-30S  gives  the  values  for  parachute  drift  for  a  wind 

of  given  velocity,  but  this  appears  to  be  based  on  the 
assumption  that  his  initial  velocity  is  zero  in  both  the 
vertical  and  horizontal  planes.  However,  if  we  are  to  deter- 
mine the  initial  navigational  error  of  the  aircraft  at  the 
instant  the  pilot  leaves  the  plane,  the  pilot  is  traveling 
with  a  velocity  equal  to  that  of  the  plane. 

To  write  an  equation  for  this  component  of  his  motion, 
we  will  consider  that  there  is  no  wind  and  the  plane  is 
flying  horizontally  at  constant  speed. 

We  will  define  the  y-axis  as  the  vertical  and  the  - 
x-axis  as  the  horizontal  in  the  direction  of  his  motion. 

Since  we  know  that  his  initial  velocity  is  that  of 
the  airplane,  and  there  is  no  force  acting  on  him  except 
that  of  air  resistance,  he  will  decelerate  to  zero  velocity. 
And  since  his  initial  velocity  is  so  great,  we  can  readily 
assume  that  the  resistance  is  proportional  to  the  square 
of  the  velocity. 
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Horizontal  Motion  Due  to  Initial  Velocity 


k  r  constant  or  factor  of  proportionality 

Vu  ss  velocity  due  to  that  of  the  airplane 

V0  -  initial  velocity 

mdVn  *  -k2V  2      . 
dt 


An 


-k  dt 
m. 


"vt" 


r- .  -*2 


v, 


m 


Jl_-  JL_»  fcfJL 


V. 


u 


Yr 


m 


1  -  St  =  vu!ai 

=   ^U  frnk2t    +    l\ 

lv0m         ; 


=  vuf 


u 


ml 


£L 


Vnk2t  +  m 


=  mV0/' 1 \ 

At     tsp,      x  -  0 


dx  = 
dt 


J  o  y 


raVo       dt 
o       V0k^t  - 


ra 


(6) 
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x  -  m  ln(Vrrk2t  +  ra}l 
k2"  Jo 


"2 

x  =  _m_Ln(V0k  t  +  m)  -  ln(m) 

k2 


x  -  m  ln/V0k  t  •»  m 
k    V    a 


X 

s     m  in 

k2 

Let 

b2 

- 

-*- 

k^ 

X 

=  b2ln 

/Vpk2t  t  l\ 


(7) 


Here  I  must  state  that  a  few  other  assumptions  have 
been  made.   These  are: 

1.  The  rotation  of  the  earth  has  been  neglected, 

2.  The  coriolis  force  has  been  neglected. 

3.  The  sphereicity  of  the  earth  has  been  neglected. 
I  feel  that  these  assumptions  would  not  materially 

affect  the  accuracy  of  the  drift.   This  is  so  because  of  the 
short  distances  involved. 

Initially,  the  formulae  were  derived  for  conditions 
giving  a  limiting  velocity  of  175  ft/sec.   This  value  was 
taken  from  Differential  Equations  with  Engineering  Applica- 
tions, page  51.  We  know  that  air  resistance  or  parasite 
drag  has  the  form  D  Z  Cd£__k    (Theory' of  Flight— von  Mises— 

oage  95).   This  corresponds  to  c  sr  Cd  r   A.  in  our  initial 

2 
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derivation.   There  we  assumed  that  the  value  of  c  was  con- 
stant.  A  quick  glance  at  this  equation  immediately  reveals 

o 

that  c  is  not-  constant  since  it  depends  on  the  air  density 

which  varies  with  altitude. 

The  airea,  A,  is  the  projected  or  frontal  area  of  the 
body  under  consideration.   In  order  to  get  the  maximum 
value  of  drift,  I  have  assumed  that  the  projected  area  of 
the  man  will  be  the  effective  area  in  all  three  cases  of 
motion.   Since  a  body  will  move  with  the  frontal  area  normal 
to  the  direction  of  movement,  we  will  use  this  in  free  fall, 
motion  due  to  the  wind,  and  motion  due  to  initial  velocity 
and  thus  have  the  maximum  drift.   This  value  will  then  be 
a  constant. 

The  density  of  the  air,  £>  ,  depends  on  the  barometric 
pressure  and  the  temperature  of  the  air.  Since  each  of  these 
vary  with  altitude,  the  density  in  turn  varies  with  altitude, 
I  will  use  the  values  of  pressure  and  temperature  for  the 
U.  S.  Standard  Atmosphere  and  plot  these  values  against  - 
altitude.   Hopefully,  I  will  be  able  to  write  an  equation 
defining  this  variation. 

The  third  factor  is  the  drag  equation  is  C^,  the  co- 
efficient of  drag.  C<2  depends  on  the  shape  of  the  body  and 
the  Reynolds  number  of  the  flow. 

rr>    ~  VL  -  Reynolds  Number 
V 

V  -  velocity 
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L  =  length 
-£/  s  '       ■  kinematic  viscosity  of  the  air 

jf   -  dynamic  viscosity  of  the  air 
/>  -  density  of  the  air 

The  dynamic  viscosity,  u  ,  varies  with  the  temperature 
but  not  with  pressure,   (Applied  Fluid  Mechanics  -  page  333) • 
Therefore  we  may  plot//  against  temperature  and  in  turn  con- 
vert it  to  altitude*  With  this  and  the  plot  of  density,  we 
may  plot  the  kinematic  viscosity,  -y  ,  versus  altitude. 

With  this,  we  may  investigate  the  variation  in  "TR 
and  thus  the  variation  in  C^. 

Data  for  plot  of  Dynamic  Viscosity  versus  Temperature, 

Ref :   Table  52  -  Handbook  of  Meteorology.  -  Berry,  Bollay, 

and  Beers  -  McGraw  Hill  -  New  York  -  1945. 

Temperature  Temperature  Dynamic  Viscosity  Dynamic  Viscosity-, 

°C         °F      gram/cm-sec.       pound-sec. /ft? 

(x  10-°)  (x  10-7) 

-70  -94  172  3*595 

-50  -53  147  3.075 

-30  -22  152  3.130 

-10         14  162  3.330 

0        32  172  3.595 

10         50  177  3.700 

20         63  131  3.730 

30        36  136  3.330 

40        104  191  3.99 
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Data  for  plot  of  Barometric  Pressure,  Temperature,  and 
dynamic  Viscosity  versus  Altitude.   Ref:  Pressure  and  Temp- 
erature -  Handbook  of  Meteorology.  Dynamic  Viscosity  -  Plot 
of   vs.  Temperature. 


Altitude  Barometric  Pressure  Temperature  Dynamic  Viscosity 


Ft.x  1000 

mm,    of  Hg 

°F 

pound- se  c^/ft? 
Cx  10~7) 

0 

760.0 

59.0 

3.735 

1 

732.9 

55.4 

3.720 

2 

706.6 

51.3 

3.707 

3 

631.1 

43.4 

3.695 

4 

656.3 

44.3 

3.630 

5 

632.3 

41.2 

3.663 

6 

609.0 

37.6 

3.643 

7 

536.4 

34.0 

3.613 

3 

564.4 

30.6 

3.579 

9 

543.2 

27.0 

3.526 

10 

522.6 

23.4 

3.473 

11 

502.6 

19.3 

3.430 

12 

433.3 

16.2 

3.396 

13 

464.5 

12.6 

3.369 

14 

446.4 

9.1 

3.343 

15 

423.3 

5.5 

3.325 

16     . 

411.3 

1.9 

3.304 

17 

395.3 

-1.7 

3.236 
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Data  continued: 

13  379.4  -5.3  3.226 

19  364.0  -3,7  3.243 

20  349.1  -12.3  3.230 

21  334.7  -15.9  3.212 

22  320.3  -19.5  3.193 

23  307.4  -23.1  3.174 

24  294.4  -26.5        "  3.157 

25  231.9  -30.1  3.140 

26  269.3  -33.7  3.125 

27  253.1  -37.3  3.111 
23  246.9  -40.9  3.100 

29  236.0  -44.5  3.090 

30  225.6  -47.9  3.033 

31  215.1  -53.5  3.0756 

32  205.3  -55.1  3.0750 

33  196.4  -53.7  3.076 

34  137.4  -62.3  3.031 

35  173.7  -65.7  3.093 

36  170.4  -67.0  3.101 

37  162.4  -67.0  3.101 
33  154.9  -67o0  3.101 

39  147*6  -67.O  3.101 

40  140.7  -67.0  3.101 

41  134.2  -67.O  3.101 
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Data  continued: 


42 

127.9 

-67.0 

3.101 

43 

122.0 

-67.0 

3.101 

44 

116.3 

-67.0 

3.101 

45 

110. 0 

-67.0 

3.101 

46 

105.7 

-67.0 

3.101 

47 

100.7 

-67.0 

3.101 

4S 

96.05 

-67.0 

3.101 

49 

91.57 

-67.0 

3.101 

50 

37.30 

-67.0 

3.101 

The  density  of  the  air  and  the  kinematic  viscosity 
may  now  be  calculated. 

where: 


Z3  r"~gfe" 


p  s  pressure 

g  -  acceleration  due  to  grav- 
ity (assumed  constant) 


R 

»  gas  constant  for  Air 

i# 

T 

-  temperature 

Altitude 

Densitv  r  / 

Kinematic  Viscositv 

Ft.  x  1000 

Slugs/ft? 
Cx  10-3) 

Ft?/sec.  x  10~4 

0 

2.3^ 

1.569 

1 

2.305 

1.614 

2 

•  2.245 

1.652 

3 

2.159 

1.711 

4 

2.115 

1.740 
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Data  continued: 

5 

2.050 

1.739 

6 

1.935 

1.S3S 

7 

1.930 

1.372 

-  * 

1.363 

1.916 

9 

1.313 

1.942 

10 

•  1.759 

1.976 

11 

1.705 

2.010 

12 

1.643 

2.061 

13 

1.595 

2.110 

14 

1.543 

2.165 

15 

1.492 

2.223 

16 

1.449 

2.234 

17 

1.400 

2.343 

IS 

1.357 

2.412 

19 

1.309 

2.435 

20 

1.264 

2.555 

21 

1.225 

2.620 

22 

1.130 

2.704 

23 

1.142 

2.765 

24 

1.103 

2.360 

25 

1.069 

2.947 

26 

1.023 

3.043 

27 

0.994 

3.130 

23 

0.960 

3.232 
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Data  continued: 

29 

0.921 

3.354 

30 

0.333 

3.475 

31 

0.364 

3.540 

32 

0.325 

3.730 

33 

0.794 

3.375 

34 

•  0.765 

4.030 

35 

0.735 

4.210 

36 

0.705 

4.405 

37 

0.671 

4.624 

33 

0.639 

4.360 

39 

0.610 

5.094 

40 

0.531 

5.342 

41 

0.555 

5.595 

42 

0.529. 

5.376 

43 

0.504 

6.165 

44 

0,431 

6.450 

^ 

0.457 

6.300 

46 

0.436 

7.121 

47 

0,416 

7.470 

43 

0.397 

7.331 

49 

0.379 

3.202 

.  50 

0.361 

3.610 

Now  that  we  have  the  various  graphs  plotted,  we  have 
to  investigate  the  variation  in  the  Reynolds  Number  to  deter* 
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mine  the  variation  in  the   drag  coefficient. 
|fc     s  Hi 

At  50,000  feet:  <p  =  S.61  x  10^ 

xj  r  0.361  x  10-3 

/«.  =.  3.101  x  10-7 

At  10,000  feet:    .  sy  z  1.976  x  10"^ 


/° 


-  1.759  x  10~3 


Zi  =  3.473  x  10~7 
The  weight  of  the  man  need  be  revised  to  account  for 

the  weight  of  the  parachute  and  the  flight  suit.  If  we 

figure  20  pounds  each  for  the  last  two  items,  mg  -  W  s  220 

pounds. 

We  must  assume  a  shape  for  the  man.  There  are  two 

possibilities: 

1.  Circular  cylinder,  d  =  1.5  feet 

2.  Ellipsoid,   L  -  6  feet,   d  -  1.5  feet 

In  both  cases,  the  Reynolds  Number  will  be  the  same. 

At  50,000  feet:     f|7  :  Vd     B     7(1.5)  =  1.733  x  103V 

1  -p  8.61x10-* 

At  10,000  feet:    Tj?    r  M    =       V(l*5)  -  7.59  x  103V 

1J         1.976  x  10"* 

When  the  man  is  wearing  a  flight  suit  with  boots  and 
helmet  and  has  a  parachute  on,  he  will  more  closely  approxi- 
mate a  circular  cyclinder.  For  this  reason,  and  the  fact 
that  there  is  very  little  data  on  ellipsoids,  I  will  use  the 
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form  of  a . circular  cylinder. 

From  WeisbergerTs  experiments,  if  we  can  keep 
TR>5  x  1CK,  we  will  have  a  constant  value  for  C^.  (Applied 
Hvdro-and-Aero-Mechanicsr  Prandtl  and  Tietjens). 

We  will  examine  IK  at  both  50,000  and  10,000  feet  to 
determine  what  velocity  we  must  have. 
At  50,000  feet:   ^  -  5  x  105  -  1.733  x  lO^V 

V  must  equal  28$   ft. /sec. 
At  10,000  feet:  TR  -  5  x  105  -  7*59  x  lO^V 

V  must  equal  65.9  ft ./sec. 

To  see  if  this  is  possible,  we  must  check  the  limit- 
ing velocity  at  each  altitude.  We  know  that  the  limiting 
velocity  is  reached  when  the  drag  equals  the  weight  since 
there  is  no  acceleration. 

With  IK  s  5  x  10^,  Cd  -  0.33  for  an  infinite  cylinder 

(Handbook  of  Engineering  Fundamentals  -  Eshbach) . 

D  -  W  3  CU/?  V2A 
2 

At  50,000  feet:   V2  -  2W   -  2(220)  (l.U) 

Cd/>A   0.33(0.361  x  10-.>H11.25) 

V2  -  47.3  x  104 

V  =  635  ft. /sec. 

We  found  that  we  must  have  a  velocity  of  2S£  ft. /sec.  so  our 

assumption  is  valid  at  50,000  feet. 

At  10,000  feet:    V2  -  2W _  2(220)  (1U.) 

C^A  ~  0.33(1.759  x  10-3)  (1125) 

V2  =  9.3  x  104 
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V  =.  314  ft. /sec. 

We  must  have  a  velocity  of  65.9  ft. /sec,  so  our 
assumption  is  also  valid  at  10,000  feet. 

This  criterion  of  TR>  5  x  10^  will  hold  true  every 
where  between  10,000  and  50,000  feet.   This  becomes  obvious 
when  examining  the  curves. 

WithTR)  5  x  10',  we  may  enter  Table  III  in  Rouse, 
page  249,  and  find  that  for  a  length-diameter  ratio)!  of  5 
(which  is  very  nearly  what  we  have),  C^  =  0,35. 

Looking  toward  the  same  problem  with  a  parachute, 
this  investigation  has  also  revealed  a  value  of  C^  =  1.33 
for  a  parachute.  (Theory  of  Flight  -  von  Mises  -  page  9$). 
This  value  is  in  agreement  with  the  values  found  in  other 
texts  for  a  hemisphere,  hollow  upstream.   The  criterion  here 
is  that  TR  >  10-*  which  is  easily  recognized. 

It  should  be  noted  that  these  values  of  critical 
Reynolds  Numbers  decrease  with  increasing  turbulence  of  the 
air.   The  values  used  herein  are  for  air  with  no  turbulence 
in  the  approaching  stream.   (Fluid  Mechanics  -  Cox  and 
Germano  -  p.242).   This  gives  us  a  factor  of  safety  since 
the  air  will  be  turbulent. 

We  have  now  covered  the  factors  with  one  exception, 
density.   Therefore, the  problem  at  hand  is  writing  an 
equation  which  describes  the  density  variation. 

From  the  plot  of  density  versus  altitude,  it  appears 
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that  we  have  a  decreasing  exponential.   This  will  give  us  an 
equation  in  the  form  of 

^   =  Ae~By 

A  plot  of  yO   vs.  altitude  on  semi-log  paper  will  give 
us  an  indication  of  whether  or  not  we  have  an  exponential 
relationship.  (See  graph  next  page). 

By  sighting  down  the  plot,  it  is  seen  that  we  do  not 
have  a  straight  line.   However,  a  straight  line  fits  the 
curve  from  0  to  32,000  feet,  and  another  from  32,000  to 
50,000  feet. 

I  would  rather  approximate  the  relationship  with  one 
smooth  curve  rather  than  make  it  discontinuous.   The  problem 
we  now  have  is  that  of  tolerances.  How  much  tolerance  can 
we  allow  and  where  should  errors  be  a  minimum? 

With  this  determined  we  only  have  to  integrate  and 
the  manTs  motion  is  defined. 
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RLA 

11  March  1964 
From  the  plot  of   vs.  y  on  semi-log  paper,  we  may 
approximate  the  curve  with  a  straight  line.   The  equation 
has  the  form 

Let    -  2.40  x  lCf5 

For  our  first  approximation,  consider 

y  ^  35,000  feet,  where   =  0.735  x  10-3 
From  this,  B  r  -3.375  x  10~5 

With  this  as  a  starting  point,  we  may  find  the  dif- 
ferences between  the  formula  value  and  the  actual  value. 
The  value  of  B  may  then  be  adjusted  to  give  differences  which 
are  acceptable. 

Several  adjustments  give  B  -  -3.44  x  10*"^.   For  this 
value,  the  following  differences  exist. 

B  =:  -3-44  x  10~5 
(Values  for  Actual  density,  formula  density  and  the  difference 
between  formula  density  and  actual  density  are  all  expressed 
as  slugs/ft.  x  10""^.  Altitude  is  in  thousands  of  feet). 

Altitude   Actual  Density   Formula  Density   Difference 


0 

2.3S 

2.40 

0.02 

5 

2.050 

2.021 

-0.029 

10 

1.759 

1.700 

-0.059 

15 

1.492 

1.433 

-0.059 

20 

1.264 

1.209 

-0.055 
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a  continued: 

2* 

1.069 

1.040 

-0.029 

30 

0.353 

O.S56 

-0.032 

35 

0.735 

0.731 

-0.004 

40 

0.551 

0.607 

0.026 

45 

0.457 

0.510 

0.053 

50 

0.361 

0.430 

0.069 

We  now  have  an  equation  for  the  density  of  air  as  a 
function  of  altitude. 

/°   =  2.W  x  HT^-M*  x  10"5y) 

Now  we  can  proceed  and  integrate  our  basic  equation 
for  free  fall.   At  this  time,  I  will  use  the  letter  h  to 
designate  altitude  and  y  will  be  used  strictly  for  vertical 
displacement.   Therefore,  y  =  h0  -  h  where  h0  is  the  alti- 
tude at  which  the  fall  began.   The  relationship  for  air  den- 
sity now  becomes 

/>  =  2.40  x  10-3e(-3'U  *  10"5h) 
Gur  primary  equation  is  now 

^   dt        ^2 

^t  CdA,  iL,  -  C.      Then     mcilm=  nig  -  CY2eBh 
rZ  dt 

Since     y  =  h0  -  h,     h  s  ho  -  y    where  y  =.  J Vdt 

mdV^  =  mg  -  CVg  eBho 
dt  e^ 
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For  a  given  altitude,  e  °  is  constant. 

Let  CeBfl°  =  D 

Since  3  is  negative,  we  have 

mdV^  =  mg  -  DVge1^ 
dt 

I  have  tried  for  -over  a  week  to  integrate  this  equa- 
tion consulting  every  differential  equation  textbook  that  I 
could  find.   This  has  not  produced  any  positive  results. 
Since  this  equation  appears  to  be  impossible  to  integrate, 
we  can  also  approximate  the  relationship  for  density  as  a 
series  of  straight  lines  of  changing  slope.  We  may  consider 
each  curve  over  a  definite  interval  and  take  three  intervals 
in  order  to  make  a  good  approximation.   This  equation  has 
the  form 

P    -  ah  +  b 
The  first  interval  we.  shall  consider  is 

(1)  0  ^  h  *  14,000  feet 
At  h  :  0,  />  :  2.37  x  10~3 

r>  bx  =  2.37  x  10~3 
At  b  a  14,000,  /o   =  l053  x  10-3 
z$>>   a2  =.  -6.00  x  10-3 

Therefore,  /O  j  =  -6.00  x  10"%  +-  2.37  x  10"3 
The  second  interval  is 

(2)  14,000  ^  h  ^  30,000  feet 
At  h  -  14,000,^  =  1.53  x  10"3 
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:=>  b2  =  1.53  x  10-3 

At  h  =  30,000,    /?  S  0.37  x  10 ^ 

'Z>   a2  =  -4.125  x  10"d 

Therefore,    y02  =  -4.125  x  10~^h  4  1.53  x  lCT3 
The  third  Interval  Is 

(3)      30,000  ^  h  ^  50,000  feet 

Ath  =  30,000,  yO  =  0.37  x  10~3 

r>  b3  3  0.37  x  10**3 

At  h  =  50,000,   p  -  0.34  x  10~3 

z$   a^  s  -2.65  x  10"S 

Therefore,    /O3  =  -2.65  x  10~%  +  0.#7  x  10~5 

With  our  density  described  in  another  form,  we  may 
again  write  the  differential  equation  for  the  manTs  motion. 

Now  we  have 

2 

mdV^  a  mg  -  Cd(ah,  +  b)VmAm 

2  2 

mdV^i  =  mg  -  CnaATnhV^  -  C^bA^Vm- 

dt  2  2 

Let  C^Am  s  Caf  then  mdVm  <*  mg  -  C'aahV^  -  CabV& 
"2  dt 

But  again,  y  =  J  tfdt  =■  h^,-  h,  and  h  -  h0  -  y 
mdVm  B  mg  -  Caa(h0  -  y)vg  -  C  bV& 

mdVnr  s  mg  -  (Ca^o  +  ca^^vm  +  CaayVJ| 

dt 
Again,  X  have  worked  several  days  on  this  equation  and 

have  been   unable  to  integrate  it. 

If  we  can't  integrate  this  equation,  or  the  one 
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involving  the  exponential  relationship  for  density,  then  we 
must  find  some  way  of  approximating  this  solution* 

Assume  for  the  present  that  we  have  a  constant  density. 
In  that  case,  our  original  integration  holds  and  we  have 

(1)  Vm -  utanhCgt/u) 

(2)  y  r   (uA)ln  coshCxt) 

A 

u2  =_g_;     k2  =  _s£;     c2  r  G&dJLmi    X   =  -&-J    X2  *  g2; 

k  m  2  u  5  " .   u. 

X2  x       g       s  k       and       X     -  k     s     c2 
~      g2/k^  u  m 


X    "     g  gk         cVlT        S      ~V  c~d/0_Am 


-Ji-  -       ^W  and     A  =  g  -j/Cd^A^ 

"      X         g^  /      2W 

Therefore,      u  -      /     2V/ 

u  K       2W 

Jl^^gC^fim 

u  2W 

In  order  for  us  to  test  the  validity  of  this  approxi- 
mation, we  must  determine  the  numerical  values  and  express 
these  constants  as  functions  of  the  density. 

u  =  ,  /  2W    s  12.7/-/P~ 

-£-  = SL =  2,535-/^ 
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,  u_,  ~        2W  =   5.02 

A  gCd/oAra  ^o 

1     =r     1    /      2W       =  0,394  : 

A  gV  cd/°Arn       V^~ 

-A.  =  gCfT/0^  r  0.199/> 
u  2W  ' 

We  may  now  rewrite  our  equations  as 

(2)     Vm  =  12^tanh(2.535V7Tt) 

(2)  y  -   5,02  In  cosh(2.535-V/T  t) 

(3)  t  =  QJt2aitcosh~1[e(o-19^oy)] 

If  we  analyze  this  method  of  solution  over  a  large 
vertical  displacement,  we  find  that  effectively,  we  are  con- 
sidering the  density  to  be  constant  in  accordance  with  our 
initial  assumption.   This  will  provide  very  inaccurate 
results.   Therefore,  if  we  are  to  assume  constant  density, 
we  must  limit  our  value  of  y  to  a  small  increment.   I  will 
arbitrarily  choose  a  value  of  1,000  feet  for  y  as  the  inter- 
val of  displacement.   This  I  feel  is  a  safe  assumption  and 
should  not  introduce  any  significant  error.   If  we  are  to 
proceed  in  this  manner,  we  must  account  for  the  initial 
velocity  at  the  beginning  of  each  of  the  intervals.   To  do 
this  we  must  return  to  our  original  derivation  where 

mdV^  =,  mg  -  c2Vm 
dt 
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At  t  s  t,  Vm  =  Vm 
At  t  -  0,  V^  =  V0 
If  we   integrate, 

tanh^CTn/ujT      s  Cg/u)t 

tanh^CVjjj/u)    z   (g/u)t  +  tanh-l(V0/u)      a^d 
Vm  s  utanh[(g/u)t  +  tanh^CVo/ufj 
Integrating  again 

7  K  Cu2/g)ln  cosh[(g/u)t  +  tanh^C^/u)] 
t  r     1    [cosh"1^7^  -  tanh^CVu)] 


t:0, 


32&rcosh-1(e0-199/,y)   -  tanh-Vvoa^ 
^"L  V     12-7 


\7 


-7 

To  pursue  this  method  of  solution,  we  must  use  the 
following  method  of  approach,  and  we  may' use  the  exact 
values  of  air  density.  We  will  consider  h0  -  20,000  feet. 
For  y  -  1,000  feet,  we  will  determine  the  time  for  free  fall 
to  19,000  feet o assuming  a  constant  air  density  equal  to  the 
value  for  19,000  feet.  We  will  also  determine  the  velocity 
at  19,000  feet  which  will  be  the  initial  velocity  which  we 
shall  use  for  the  interval  between  19,000  and  13,000  feet. 
Results  of  these  calculations  follow. 


h 

y° 

x  10~3 

m. 

Vo 

t 

19,000 

1.309 

223 

0 

3.16 

13,000 

1.357 

224 

223 

4.43 

17,000 

1.400 

222 

224 

4.50 

16,000 

1.449 

221.5 

222 

4.51 
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In  this  case,  t  does  not  give  the  time  of  free  fall 
except  for  the  first  1,000  feet.   It  appears  that  here  the 
man  has  reached  the  limiting  velocity.  Since  the  velocity 
changes  very  little,  we  may  divide  the  average  velocity  for 
the  interval  into  1,000  feet  to  find  the  time  of  fall  for 
that  interval. 

If  we  started  the  same  procedure  at  say  h0  -  25,000 
feet,  and  at  17,000  feet  had  a  velocity  of  222  ft. /sec,  we 
could  be  reasonably  sure  that  our  approximation  is  correct. 
Let  h0  -  25,000  feet. 


h 

P 

x  1CT-* 

\ 

Vo 

t 

24,000 

1.103 

227 

0 

3.14 

23,000 

1,142 

22^ 

227 

4.39 

22,000 

1.1S0 

226.5 

22^ 

4.41 

The  tendency  of  these  calculations  is  to  verify  each 
other.  Our  first  thought  would  then  be  that  this  approxi- 
mation is  correct.  But,  let  us  look  ait  the  limiting  veloci- 
ties at  the  different  altitudes. 

h  -  17,000  feet 

Y2  -   2¥ 
~   cd/>Am 

v"  -  339  ft./sec. 
h  -  24,000  feet,  V  -  3#2  ft./sec. 
By-  our  equation  in  the  approximation,  consider  the 
.interval  of  fall  from  50,000  to  49,000  feet. 
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49,000  ^  h  ^  50,000;    V  -  245  ft. /sec. 
4S,000  *  h  *  49,000;    ?  =  222  fto/sec* 

The  limiting  velocity  at  43,000  feet  is: 

V  =  636  ft. /sec. 

In  all  three  cases,  we  notice  quite  a  disparity  be- 
tween the  velocity  as  determined  by  our  equations  and  the 
limiting  velocity.   From  this,  I  think  we  must  accept  the 
approximations  as  being  invalid. 

In  order  to  determine  the  time  of  free  fall,  another 
method  of  approximation  could  be  used.  We  could  assume  that 
the  man  has  reached  his  limiting  velocity  and  plot  these 
values  of  limiting  velocity  versus  altitude.   By  taking  the 
average  velocity  over  a  1,000  foot  interval,,  we  could  find 
the  time  taken  to  traverse  that  interval.  Now,  suppose  that 
this  would  give  us  results  without  appreciable  error.  We 
would  then  have  no  means  to  determine  his  horizontal  motion 
because  we  have  no  knowledge  of  whether  or  not  he  will  reach 
his  limiting  velocity  in  the  horizontal  direction. 

For  an  accurate  solution  then,  we  must  integrate  our 
equations  of  motion.   I  v/ill  now  write  these  equations  and 
define  all  of  the  terms  in  summary  form. 

Our  first  set  of  equations,  (l),  (2),  and  (3),  are  for 
the  condition  of  free  fall.   Bear  in  mind  that  the  equation 
for  the  density  of  air  may  be  either  the  exponential  rela- 
tionship or  the  linear  relationship. 
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Cl)  Motion  vertically  downward. 

f  eBl\) 
mdV^  -  mg  -  CW^r   (v^Am 
dt         2  (ah  +  bj 

C2)  Horizontal  motion  due  to  wind. 


2  £ah  4  bj 


(3)   Horizontal  motion  due  to  initial  velocity. 


fln 


mdY^  =  -_CdV  ~or  \V^Am 
dt      2  /ah 


or  \Y-l 
f  b) 

It  is  easily  seen  that  the  results  of  equation  (l) 
relating  the  displacement  in  free  fall  must  be  an  input  into 
equations  (2)  and  (3)  in  the  relationship  for  density. 

The  various  terms  are  defined  below. 

m  -  mass  of  the  man  -  W  in  slugs 

S 

W  -  weight  of  the  man  in  pounds 

g  s  acceleration  due  to  gravity  s.  32.2  ft./sec. 

Vm  -  velocity  of  the  man  in  the  vertical  direction  in 

ft. /sec. 
V^  -  velocity  of  the  man  in  the  horizontal  direction 

due  to  the  wind  in  ft. /sec. 
^u  s  velocity  of  the  man  in  the  horizontal  direction 

due  to  the  initial  velocity  of  the  aircraft  in 

ft. /sec. 
cd  s  coefficient  of  drag  =0.35 
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Am  -  projected  area  of  the  man  r  1125  square  inches 
\r  =  velocity  of  the  wind  in  ft. /sec.        «. 
v  -  weight  density  of  the  air  r  p  g  r  g-f       or 

'         [3&    +   b 

^Oo  :  2.40  x  ICT^slugs/ft? 

B  r  «3.44  x  10~5  ftT1 

h  -  altitude  in  feet 

h0  -  altitude  of  bail-out 

y  -   vertical  displacement  in  feet 

h  =  h0  -  y 

a  and  b  are  defined  for  the  following  intervals: 

0  £  h  ^  14,000     a  -  -6,00  x  10~f 

b  *  2.37  x  10~-* 

14,000  ^  h  ^  30,000     a  =  -4,125  x  107^ 

b  s  1.53  x  10~3 

30,000  *  h  ^  50,000     a  -  -2.65  x  10"^ 

b  s  0.37  x  10~? 

We  must  now  look  at  the  equations  for  fall  with  a 
parachute.  Basically,  there  is  no  difference  except  for  the 
values  of  0^  and  the  projected  area. 

For  motion  in  the  vertical  direction,  we  may  use  a 
value  of  Cd  -  1.33  (Theory  of  Flight  -  von  Mises  -  page  93} 
Tor  a  parachute.   This  value  is  in  agreement  with  those 
found  in  other  textbooks  for  a  hemisphere,  hollow  upstream, 
for  a  Reynolds  number  greater  than  10^.   Hereafter,  I  will 
refer  to  this  drag  coefficient  with  the  symbol  C^  . 

Since  we  have  assumed  the  shape  of  a  hemisphere,  the 
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projected  area  for  motion  in  the  vertical  direction  will  be 
the  area  of  a  circle  of  a  radius  equal  to  fourteen  feet. 
Therefore,  Ay  -TTCl4)2  s  615.4  ft? 

For  motion  in  the  horizontal  direction,  the  projected 
area,  Ah  -  |<rr(l4)2  -  307*7  ft.2  This  is  the  area  of  a  semi- 
circle of  a  radius  of  fourteen  feet. 

The  determination  of  the  drag  coefficient  for  side- 
wise  motion  presents  a  problem  since  there  have  been  no 
experiments  performed  for  a  hemisphere  with  its  major  axis 
perpendicular  to  the  direction  of  flow.   In  the  vertical 
direction,  C^  -  1.33  with  the  hollow  end  upstream.  For  a 
hemisphere  with  its  hollow  end  downstream,  C^  =  0.34.  After 
checking  in  many  fluid  mechanics  textbooks  and  studying  the 
various  drag  coefficients  for  different  bodies  under  several 
conditions  of  flow,  and  studying  theoretical  flow  patterns 
and  pictures  of  actual  flow  patterns,  I  will  assume  a  drag 
coefficient  for  the  parachute  in  the  horizontal  direction. 
The  value  which  I  have  chosen  is  C^  -  0,6*0. 

Now,  we  may  write  our  equations  of  motion  for  para- 
chute fall. 

(4)  Vertical  motion. 


l^-^U0^^ 


(5)  Horizontal  motion  due  to  wind. 

md7h   -    S  Yhh  -  £dfJ /°°or     (  V2Ah 
dfc  2g  2    [ah+  bj 
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(6)   Horizontal  motion  due  to  initial  velocity. 

mdV„   -  -cdh  i^or     Vvf.Ah 
dt  2    (ah  +  bj 

C^y      B     1*33 

Ay  z  615.4  ft? 

cdh  -  o.ao 

Ah  *  307.7  ft? 

The  rest  of  the  terms  are  the  same  as  defined  for 
equations  (l),  (2),  and  (3). 

In  all  of  the  above  equations,  all  must  be  integrated 
twice  in  order  to  give  a  displacement  in  the  particular 
direction.  Equations  (l)  and  (4)  must  be  solved  for  t  by 
first  integrating  to  find  the  displacement  as  a  function  of 
time  and  then  transposing  to  find  time  as  a  function  of  dis- 
placement.  If  we  call  the  time  of  free  fall  given  by  equa- 
tion (l)  t^,  then  t^  substituted  into  equations  (2)  and  C3) 
gives  the  displacements  in  the  horizontal  direction.   In  the 
integration  of  equations  (5)  and  (6) ,  consideration  must  be 
given  to  the  fact  that  there  is  an  initial  velocity.  This 
initial  velocity  is  the  value  given  by  the  first  integration 
(for  V)   of  equations  (2)  and  (3)  at  time  t-^* 

By  summing  up  all  of  the  horizontal  values  considering 
the  direction  as  well  as  the  magnitude  given  by  our  equations, 
over  the  period  of  time  from  t0  to  tj  (free  fall)  and  from 
tj  to  t^   (parachute  fall),  we  have  the  parachute  drift. 
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I  would  now  like  to  write  this  problem  in  the  form  of 
a  program  flow  chart.  When  the  six  equations  of  motion  are 
integrated,  these  may  be  written  as  additional  sub-routines 
and  added  to  this  program  with  a  predefined  connector,  I 
will  also  omit  all  machine  instructions  since  all  processes 
are  mathematical  operations, 

A.  few  new  terms  were  added  in  the  flow  chart  which 
should  be  defined. 

#-  the  direction  in  degrees  toward  which  the  wind  is 
blowing. 

0    -  the-  course  of  the  aircraft  when  the  man  bails  out. 

Q   -  the  resultant  direction  of  parachute  drift. 
After  all  of  this  work,  I  feel  that  I  should  explain 
why  I  have  put  this  much  time  in  what  is  seemingly  a  minor 
matter.  In  a  recent  study  of  the  feasibility  of  the  use  of 
a  personnel  locator  beacon,  the  following  facts  were  brought 
to  light  in  a  study  of  Naval  aircraft  accidents.  The  study 
showed  that  over  9$/£  of  all  incidents  occurred  within  60. 
miles  of  a  rescue  facility.  This  is  within  the  range  of  a 
good  DF  fix.  Assuming  that  we  know  the  location  of  the  air- 
craft when  the  person  leaves  his  craft,  we  should  be  able  to 
find  him  as  quickly  as  possible.   By  accounting  for  the 
variables  discussed  in  the  analysis  of  the  problem,  we  may 
determine  where  he  will  be  in  the  water  with  reasonable 
accuracy.  In  cases  where  the  sea  and  air  temperature  are  . 
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critical,  or  where  the  man  is  injured,  it  is  absolutely 
necessary  that  we  rescue  him  as  quickly  as  possible.   It  is 
apparent  that  this  same  solution  is  applicable  over  land 
areas. 

This  analysis  of  his  motion  in  the  air  will  make  this 
possible  upon  the  integration  of  his  six  equations  of  motion. 
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APPENDIX    D. 

EXTRACTS  FROM  THE  STUDY: 
"FREE  FALL  AND  PARACHUTE  DRIFT" 

"Free  Fall  and  Parachute  Drift"  was  an  undergraduate  study  performed 
by  Cadet  First  Class  Darvy  M.  Cohan  at  the  U.  So  Coast  Guard  Academy, 
New  London,  Connecticut,   in  the  spring  of  1965.    Cohan  expanded  and 
refined  the  work  done  by  Armacost  and  Saunders  which  appears  in  the 
preceding  Appendix.     Cohan  expanded  the  problem  to  include  the  effect 
of  winds  at  three  different  altitudes  blowing  with  different  forces  and  in 
different  directions.     He  also  considered  the  possibilities  of  the 
parachute  opening  at  three  different  altitudes:  (1)  the  altitude  at  which 
the  aircraft  was  flying,   (2)  automatically  at  14,  000  feet,  and  (3)  at  the 
surface  of  the  Earth  (parachute  fails  to  open). 

Cohan  developed  three  computer  programs,  one  mainline  program 
and  two  subroutines,  all  of  which  are  included  in  this  Appendix.    A  brief 
description  of  each  program  precedes  the  program  itself.     The  flow 
diagrams  are  taken  directly  from  the  original  paper. 

VECTOR  SUM  PROGRAM 
Cohan's  entire  program  hinged  upon  the  development  of  a  program 
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to  sum  vectors  in  three  dimensions.     The  path  of  motion  represented  a 
resultant  vector  which  made  an  angle  9  with  the  horizontal.     Using  sine 
and  cosine  functions  of  that  angle  he  determined  the  vertical  and 
horizontal  distances  a  man  would  travel  in  a  specified  time  interval, 
given  the  velocity  component,   V,  within  the  air  mass,     The  total 
horizontal  displacement  relative  to  the  air  mass  is  the  sum  of  the 
distances  for  all  time  intervals  and  is  represented  by  DTR.     By  applying 
the  horizontal  motion  of  the  air  mass  the  final  resultant  motion  with 
respect  to  the  Earth's  surface  is  obtained. 

The  theory  behind  his  vector  sum  program,   VSUM,   was  this: 

1.  Establish  one  vector  as  the  major  axis  used  by  the  computer. 

2.  Build  a  triangle  using  the  angle  between  two  vectors.     One  leg 
would  then  be  the  first  vector  plus  the  projection  of  the  second  onto 
the  major  axis.     The  other  leg  would  be  the  projection  of  the  second 
vector  upon  the  orthogonal  axis.     Dividing  the  first  leg  into  the  second 
leg  produces  the  tangent  of  the  angle  of  the  resultant  vector. 

3.  Design  a  dummy  system  to  work  with  the  azimuth  angles  found  in 
the  real  world.  Convert  the  azimuth  angles  to  the  dummy  system,  solve 
the  problem,  and  convert  back  to  azimuth  angles. 
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SUBROUTINE  ORDER 

In  considering  the  three  possible  parachute-opening  altitudes, 
Cohan  ended  up  with  three  possible  vectors  at  the  Earth's  surface,   DISA, 
DISB,   and  DISC  with  their  associated  angles,  ANGA,  ANGB,   and  ANGC. 
His  SUBROUTINE  ORDER  takes  the  three  vectors,  arranges  them  from 
maximum  to  minimum  and  then  produces  the  vector  average  of  the  three. 

MAINLINE  PROGRAM 

The  mainline  program,   called  SUBROUTINE  PARAD,   uses  the  two 
preceding  subroutines  and  the  basic  equations  found  in  the  study  by 
Armacost  and  Saunders.     The  Euler  method  of  statistical  integration  was 
applied  for  time  intervals  of  two  seconds  in  the  free  fall  phase  and  four 
seconds  in  the  parachute  drift  phase.     The  exponential  form  of  the 
equations  involving  air  density  were  employed. 
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ABSTRACT 

The  available  methods  for  estimation  of  wind  currents,  mass 
transport  velocity  by  waves  and  permanent  flow  (thermohaline 
gradient  current)  are  briefly  summerized  and  a  simplified  computer 
approach  is  outlined. 

The  computed  synoptic  surface  currents  are  compared  with  monthly 
mean  current  charts  and  with  surface  wind  conditions.    This  analysis 
indicates  that  the  surface  currents  are  greatly  wind-driven.    A  detailed 
verification  procedure  which  will  use  the  observed  changes  in  sea 
surface  temperature  is  outlined. 

The  use  of  the  synoptic  current  fields  for  computation  of  divergence 
and  convergence  and  the  resulting  changes  in  subsurface  thermal 
structure  is  described.    The  relative  importance  of  the  synoptic  surface 
currents  in  ship  routing,  rescue  operations  and  other  practices  in  reviewed. 
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1.      INTRODUCTION 

A  number  of  naval,  fisheries  and  other  martime  operations  require 
a  knowledge  of  the  direction  and  speed  of  surface  currents,  as  well 
as  their  past  and  near  future  behavior.    The  Fleet  Numerical  Weather 
Facility  (FNWF)  at  Monterey,  California,  became  interested  in  ocean 
currents  primarily  because  of  their  importance  in  Anti -Submarine  (ASW) 
applications. 

Large  variations  have  been  observed  in  thermocline  depth  which 
cannot  be  explained  by  mechanical  or  convective  mixing.    These  changes 
exhibit  cycles  which  correspond  closely  to  the  evolution  of  synoptic 
weather  patterns  over  the  ocean*     It  is  quite  clear  that  current  atlases 
cannot  be  used  to  predict  thermocline  depth  when  considerable  change 
can  occur  in  a  period  of  a  few  days.    What  is  needed  for  this  particular 
problem  are  daily  current  analyses  and  prognoses. 

Navigators  also  undoubtedly  find  that  atlases,   monthly  mean  charts 
and  the  like  frequently  do  not  give  an  exact  enough  answer  to  the  question: 
what  is  the  current  at  a  given  point  in  space  and  time?    The  purpose  of 
this  paper  is  to  report  on  an  attempt  to  compute  surface  current  flow 
on  a  quasi-synoptic  schedule  and  to  show  some  preliminary  results. 
If  these  results  appear  to  be  of  use  to  navigation,  means  will  be  found 
to  accomplish  dissemination. 
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2.      BACKGROUND 

Before  selecting  and  justifying  an  approach  which  would  be  simple 
but  still  give     some  hope  of  yielding  useful  results,  it  was  necessary 
to  screen  and  evaluate  a  voluminous  amount  of  literature  on  the  subject. 
Fortunately,  a  great  part  of  this  review  had  been  done  recently  by 
Laevastu  (1962). 

Most  theoretical  approaches  have  been  mainly  concerned  with 
explaining  the  general,  more  permanent  features  of  the  horizontal 
circulation  patterns  (see,   e.g.,   Robinson  1963).     The  Russians  have 
recently  applied  correlation  theory  in  an  attempt  to  forecast  detailed 
current  changes  from  a  known  field;  however,  our  knowledge  of  the 
initial  state  (and  particularly  its  derivatives)  is  often  rather  sketchy. 
Actually,  the  ocean  responds  quite  rapidly  to  hourly  and  daily  changes 
in  driving  forces,  and  currents  are  known  to  be  variable  in  space  and 
time  (Knauss  1960), 

These  considerations  dictated  use  of  a  method  which  would  account 
for  fairly  rapid  response  and  would  stand  up  to  daily  verification.     Many 
attempts  at  current  prediction  have  been  disappointments  because  of 
oversimplifications  resulting  from  the  assumptions  made.     It  seems 
logical  to  separate  the  total  current  into  its  elementary  components  to 
see  which  should  be  neglected,  which  can  be  simplified,   etc.    This  is 
the  attack  which  has  been  followed  in  this  investigation, 
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3.      COMPONENTS  OF  SURFACE  CURRENTS 

Surface  currents  are  caused  and  influenced  simultaneously  by  a 
number  of  forces  which  vary  independently  from  each  other  in  space 
and  time.    If  one  neglects  the  special  effects  due  to  variation  in 
depth,  coastal  configuration,  runoff,  etc. ,  the  current  vector  at  a 
given  location,  time  and  depth  below  the  surface  (^XVZv )  can  be 
given  as  the  resultant  of  the  following  components: 

wxyzt  =  wc  +  ww  +  wi  +  wt  (1) 

where  Wc  is  the  permanent  flow  (thermo-haline  gradient  current  or 
"characteristic  current"  as  used  by  Palmen  (1930)  and  Hela  (1952)), 
Ww  is  the  current  due  to  transport  by  wind  and  waves,  W^  is  the 
periodic  part  of  the  inertia  current  and  Wj.  is  the  periodic  part  of 
the  tidal  current. 

The  computations  which  will  be  described  here  cover  a  period  of 
24  hours,  and  it  will  be  assumed  that  semidiurnal  and  diurnal  tidal 
components  will  equal  out  and  can  be  neglected.    In  addition,  inertial 
eddies  will  not  be  considered  because  the  available  quantitative  infor- 
mation about  their  behavior  does  not  warrant  their  inclusion  in  this 
simple  technique.    The  two  components  which  this  study  will  attempt 
to  evaluate  are  thus  the  "characteristic"  or  permanent  transport  and 
the  transport  due  to  wind  and  waves. 
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4.      THE  CHARACTERISTIC  COMPONENT 

The  characteristic  c/brnponent  is  directly  related  to  density 
gradients  caused  by  areal  differences  in  heating -cooling  and 
evaporation-precipitation.    Although  what  we  usually  call  "the 
permanent  flow"  is  strongly  controlled  by  the  large  -scale,   more- 
or-less  stationary  wind  systems,  only  the  thermo-haline  influences 
are  included  in  that  component  here.    Wind  and  wave  effects  will 
be  lumped  into  one  computation  to  be  discussed  later. 

Several  workers  have  found  (e.g.  Yasui  1957)  that  there  is  a 

close  correlation  between  ocean  temperature  distribution  and 

dynamic  depth  anomalies.    Neglecting  salinity,  one  can  apply 

the  well-known  meteorological  thermal  wind  relationship  in  the 

ocean  if  one  knows  the  mean  temperature  of  the  layer  between  the 

^surface  and  some  level  of  zero  current  velocity.    The  characteristic 

current  is  then  given  by: 

Wc    =      ?*f       vT    x    |K  (2) 

f  f 

where  g  is  gravitational  acceleration,  f  is  the  Coriolis  parameter, 
T  is  the  mean  temperature  above  the  level  of  zero  current,    ya 
is  the  depth  to  zero  current,  and  fK  is  the  unit  vertical  sector. 
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The  determination  of  representative  mean  temperature  (  T  ) 
is,  of  course,  the  critical  factor  in  this  part  of  the  problem.    The 
temperature  structure  of  the  ocean  is  certainly  not  constant,   partic- 
ularly closer  to  the  surface;  therefore,   semi-synoptic  temperature 
fields  should  be  used  if  possible,,    The  only  place  where  sufficient 
data  is  available  for  reliable  analysis  on  a  daily  basis  is  at  the 
surface  so  it  was  decided  to  use  the  FNWF  Sea  Surface  Temperature 
(SST)  analyses  based  on  84  hours  of  reported  ship  engine  injection 
temperatures  at  the  top  of  the  layer.    In  order  to  include  a  part  of 
the  deeper  temperature  structure,  the  SST  field  is  presently  combined 
with  a  climatological  field  at  200  meters  depth  to  obtain 

*   -    Kl  Tsfc    +    K2T200  (3) 

Finally,  this  field  is  modified  empirically  in  areas  where  salinity 

considerations  are  known  to  be  important  (Cyashio,  Greenland, 

Labrador  currents).    This  in  effect  corrects  the  ocean  temperatures 

for  salinity  much    as    the  meterorologist  corrects  atmospheric  temperatures 

for  moisture  content  when  he  uses  the  concept  of  "virtual  temperature.  " 
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5.      THE  WIND  COMPONENT 

According  to  Ekman  (1906),  the  direction  of  the  wind  current 
at  the  surface  is  45°  to  the  right  of  the  wind  in  the  Northern 
Hemisphere  and  this  angle  increases  with  depth.    Recent  inves- 
tigations reveal  that  the  deflection  is  more  nearly  12=20  degrees, 
being  larger  and  more  irregular  at  the  lower  wind  speeds  (possibly 
because  of  the  increased  importance  of  other  components)  and 
smaller  at  higher  wind  speeds.    As  the  surface  wind  is  about  the 
same  angle  to  the  left  of  the  geostrophic  wind,  it  is  assumed  herein 
that  the  direction  of  the  wind  current  is  the  direction  of  the  geostrophic 
wind. 

Numerous  empirical  studies  have  indicated  use  of  a  single  factor 
to  relate  surface  current  speed  to  wind  speed.    The  formula  of 
Witting  (1909)  appears  to  agree  well  with  available  data  and  further 
allows  approximate  incorporation  of  mass  transport  by  waves  in  a 
simple  expression; 


W      =   K0 

w  3    V       g 


-yWg"  (4) 

where  Wg  is  the  mean  geostrophic  wind  speed  for  a  24-hour  period. 
In  the  present  work  it  is  assumed  that  the  current  is  relatively 
uniform  and  unidirectional  in  the  turbulent  mixed  layer  down  to  the 
thermocline  (or  about  200  meters).     The  mass  transport  of  the  waves, 
however,  would  modify  this  picture  as  it  decreases  exponentially 
with  depth  (Masch  1962).    Therefore,  if  W    is  In  meters/sec  and 

182 


Ww  in  cm/sec,  K^  is  taken  to  be  3.  3  for  surface  currents  (ship  routing 
and  drift  computations)  and  2.  2  for  the  average  current  down  to  the 
therraocline  (convergence/divergence  computations).    Obviously, 
there  is  a  time  lag  between  the  change  of  the  wind  and  response  of 
the  sea.    This  lag  seems  to  be  shorter  than  previously  believed, 
however,  and  is  partially  minimized  by  the  24-hour  averaging. 

Since  all  computations  are  carried  out  in  the  standard  FNWF 
grid  system,  u  and  v  current  components  are  determined  at  approx- 
imately 200  nautical  mile  intervals  for  all  Northern  Hemsphera, 
ocean  areas.    From  these  components,  direction  and  total  transport 
(nm/day)  fields  are  determined  and  stored  on  magnetic  tape  for  later 
output  in  chart  form  or  as  special  messages  giving  the  currents  at 
specified  latitude/longitude  intersections. 
6„      RESULTS 

Figure  1  is  a  hand  analysis  of  one  of  the  first  current  computations 
made  on  a  synoptic  schedule  (18  GMT  16  November  1964).    The  con- 
tours represent  total  current  transport  in  nautical  miles  per  day; 
direction  arrows  have  been  plotted  in  the  most  significant  current 
systems.    One  can  clearly  distinguish  such  well-known  features  as  the 
Gulf  Stream,  Sargasso  Sea,  Labrador  Current,  Kuroshio  and  Cyashio. 
The  low-latitude,  westerly  return  flow  which  results  primarily  from 
the  "wind  component"  term  is  well  defined  in  both  the  Atlantic  and 
Pacific.    A  narrow  equatorial  countercurrent  was  obtained  as  a  result 
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of  the  200~meter  temperature  structure  used  in  the  characteristic 
component. 

Figures  2a  and  b  are  synoptic  current  charts  for  the  Pacific  and 
Atlantic,  respectively,  which  have  been  drawn  automatically  on  an 
incremental  x-y  curveplotter0     Each  chart  requires  approximately 
one  minute  to  complete  and  is  of  sufficient  quality  that  it  can  be 
used  immediately  for  radio-facsimile  transmission. 

Figure  3  is  a  climatological  current  chart  for  winter.    It  can  be 
seen  that  many  of  the  most  impprtant  features  are  correctly  depicted 
in  this  approach  both  in  location  and  intensity. 

The  problem  of  automatic  plotting  of  direction  arrows  on  these 
charts  has  not  yet  been  solved.     However,  a  possible  substitute 
has  been  found  and  is  now  being  programmed  for  numerical  testing. 
Since  u  and  v  current  components  are  available  in  field  form  at  all 
grid  points,  it  is  believed  a  stream  function  field  can  be  determined 
by  a  relaxation  solution  Of  the  Poisson  equation 

0             a  v         5u 
Vzy>    . (5) 

d  x  9  y 

This  would  permit  plotting  of  a  second  set  of  lines  (^0  which 
would  everywhere  parallel  the  direction  of  flow. 
7.      VERIFICATION 

A  synoptic  current  chart  would  be  of  little  value  if  it  could  not 
be  verified  and  the  computational  scheme  tuned  as  required.    Direct 
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current  measurements  in  the  open  ocean  are  too  few  and  drift 
calculations  made  from  navigational  fixes  are  frequently  inaccu- 
rate in  weak,  current  areas,   so  it  is  difficult  to  make  a  direct 
evaluation.    It  has  been  necessary,  therefore,  to  resort  to  indirect 
means  which  are  susceptible  to  verification  on  a  synoptic  basis. 

Sea  surface  temperature  (SST)  is  the  only  oceanographic  element 
which  permits  a  reasonable  complete  synoptic  analysis  on  a  hemi- 
spheric scale.    Such  analyses  are  made  twice  daily  at  FNWF 
Monterey  (Wolff  1964),  and  their  resolution  is  such  that  SST  temper- 
ature changes  can  be  determined  for  periods  of  24,   48  hours,  etc. 
From  these  changes  will  be  subtracted  the  local  changes  computed 
from  air/sea  heat  exchange  equations.    If  the  remainder  correlates 
well  with  the  advective  change  indicated  by  W       t  •  V  SST,  the 
computed  currents  can  be  assumed  to  be  reasonably  correct. 

This  method  of  verification  is  now  being  programmed  and 
numerical  results  are  not  yet  available.    Subjective  study  of 
SST  change  charts  and  corresponding  current  charts  does,  however, 
indicate  that  the  approach  described  here  is  useful.    It  is  evident 
that  the  wind  component  term  predominates  in  many  areas,  and 
that  it  is  this  term  which  is  mainly  responsible  for  the  rapid 
response  of  sea  surface  temperature  changes  in  the  ocean. 

There  are  a  number  of  modifications  which  must  undoubtedly 
be  made  to  this  program;  it  is  hoped  these  will  be  uncovered 
during  the  verification  period.    One  obvious  question  is  -  what 
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effect  does  thermocline  depth  itself  have  upon  the  surface  current 
speeds  above  the  thermocline? 
8.      APPLICATIONS 

The  surface  current  program  was  initiated  primarily  to  determine 
divergence  and  convergence  and  the  accompanying  up  and  down 
movement  of  thermocline  depth.    Furthermore,  the  results  will  be 
used  for  forecasting  the  advective  part  of  sea  surface  temperature 
changes. 

Over  large  parts  of  the  oceans  the  currents  have  little  direct 
effect  on  navigation.    In  some  areas,  however,  they  should  be 
taken  into  consideration  in  Optimum  Ship  Routing.    Charts  of  this 
type  should  also  prove  useful  in  the  prediction  of  ice  movement 
and  in  any  rescue    operations. 

It  is  planned  to  make  these  oomputations  on  a  daily  synoptic 
schedule  (probably  to  06  and  18  GMT),   and  they  could  be  trans- 
mitted from  Fleet  Weather  Centrals  in  either  facsimile  or  special 
message  format  if  such  is  desired.    Groups  such  as  the  Institute 
of  Navigation  may  determine  that  there  are  applications  in  navigation 
which  could  make  use  of  these  products. 
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APPENDIX    F. 

SURFACE  DRIFT  ELLIPSE  PROGRAM 

The  Surface  Drift  Ellipse  Program  which  comprises  this  Appendix 
was  developed  by  Donald  A.   Burns  of  the  U.  S.  Naval  Oceanographic 
Office.    The  specific  instructions  for  using  this  program  may  be 
obtained  by  writing  to  Mr.   Burns,  Code  3800,   U.  S.  Naval  Oceanographic 
Office,  Washington,   D.  C#/   20390o     The  basic  logic  of  this  type  of 
program  may  be  found  in  A  Statistical  Rose  Program,  SP-64,  available 
from  the  Evaluation  Branch,  Oceanographic  Analysis  Division,   Marine 
Sciences  Department,  U.  S.  Naval  Oceanographic  Office,  Washington, 
D.  C. 
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